FREQUENCY  CONTROL  OF  A  C02 
LASER  USING  STARK  CELL  STABILIZATION 


Daniel  Morris  Brown 


Thesis 

380993 


NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


. 

THESIS 

FREQUENCY  CONTROL  OF  A  C02 

LASER  USING  STARK  CELL  STABILIZAT 

ION 

by 

Daniel  Morris  Brown 

June  1979 

The 

sis  Advisor:                 A.  W. 

Cooper 

Approved  for  public  release;  distribution  unlimited 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (**tt  Data  Entorod) 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


1.    report  NUMICH 


2.  GOVT  ACCESSION  NO 


3.     RECIPIENT'S  CATALOG   NUMBER 


4.     TITLE  (ond  Submit) 

FREQUENCY  CONTROL  OF  A  C02  LASER 
USING  STARK  CELL  STABILIZATION 


S.     TYPE  OF   REPORT  *  PERIOO  COVERED 

Master's  Thesis; 
June  1979 


»   PERFORMING  ORG.  REPORT  NUMBER 


7.  AuTHORfaJ 

Daniel  Morris  Brown 


•  .  CONTRACT  OR  GRANT  NUMBERft.) 


t.  PERFORMING  ORGANIZATION  NAME  AND  AOORESS 

Naval  Postgraduate  School 
Monterey,  California   93940 


10.     PROGRAM   ELEMENT.  PROJECT,    TASK 
AREA  *   WORK  UNIT  NUMBERS 


II.     CONTROLLING  OFFICE  NAME  ANO  AOORESS 

Naval  Postgraduate  School 
Monterey,  California   93940 


12.     REPORT  OATE 

June    19  79 


IS.     NUMBER  OF  PAGES 


121 


14.     MONITORING  AGENCY  NAME  *   AOORESS<(/  ditiotont  from  Controlling  OtUeo) 


Naval  Postgraduate  School 
Monterey,  California   93940 


IS.     SECURITY  CLASS,  (ot  thf  rinort) 

Unclassified 


<■«.     OICLASSIFICATIOn/OOWNGRAOING 
SCHEDULE 


16.     DISTRIBUTION   STATEMENT  (of  thlm  Komort) 


Approved  for  public  release;  distribution  unlimited. 


17.     DISTRIBUTION  STATEMENT  (ot  tho  amotrmct  ontmrod  In  Block  20,  It  dltloront  feaai  Report) 


1».     SUPPLEMENTARY  NOTES 


1t.     KEY  WORDS  (Contlnuo  on  rororoo  *tdo  It  nocooonry  fd  Identity  by  block  nimnbor) 

Laser 

Stabilization 

Stark 


20.     ABSTRACT  (Contlnud  on  rovoroo  tido  It  nocoa—ry  and  Identity  by  alock  mmmkor) 

Preliminary  studies  have  been  carried  out  for  the  frequency 
and  intensity  stabilization  of  the  output  of  a  Sylvania  Model 
941  Carbon  Dioxide  laser,  by  feedback  control  of  resonator 
length.   Fluctuations  due  to  discharge  current  and  temperature 
and  flow  rate  of  coolant  were  studied.   Current  regulation 
of  the  power  supply  was  found  to  be  necessary  for  intensity 
stability.   Cavity  length  tuning  allowed  vibration-rotation 


DD  ,  KE.  1473 

(Page  1) 


EDITION  OP  1  NOV  «•  IS  OBSOLETE 

S/N    0102-014-4601   | 

1 


SBCURITY  CLASSIFICATION  OF  THIS  PAOE  (Whon  Dmta  gntorod) 


line  selection,  with  a  strong  hysteresis  effect.   A  computer 
analysis  of  error  signal  production  for  a  pyroelectric 
detector  and  for  a  HgCdTe  detector  indicated  that  phase 
sensitive  detection  of  the  first  harmonic  signal  component 
gave  the  typical  frequency  discriminant.   Operation  of  a 
servo  loop  integrator  for  the  feedback  loop  has  been  con- 
sidered.  A  suitable  Stark  Effect  absorption  cell  for  feed- 
back signal  generation  has  been  constructed,  but  has  not 
been  applied  to  the  laser. 


DD  Form   1473 

.  1  Jan  73  

S/N  0102-014-6601  2  mcu«mtv  claudication  o'  ▼*•«  *»»otr»*««  o* 


Approved  for  public  release;  distribution  unlimited 

FREQUENCY  CONTROL  OF  A  CO 
LASER  USING  STARK  CELL  STABILIZATION 

by 


Daniel  Morris  Brown 

Lieutenant,  United  States  Coast  Guard 

B.S.,  United  States  Coast  Guard  Academy,  1975 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 


MASTER  OF  SCIENCE  IN  PHYSICS 


from  the 

NAVAL  POSTGRADUATE  SCHOOL 
June  1979 


ABSTRACT 

Preliminary  studies  have  been  carried  out  for  the 
frequency  and  intensity  stabilization  of  the  output  of  a 
Sylvania  Model  941  Carbon  Dioxide  laser,  by  feedback 
control  of  resonator  length.   Fluctuations  due  to  discharge 
current  and  temperature  and  flow  rate  of  coolant  were 
studied.   Current  regulation  of  the  power  supply  was  found 
to  be  necessary  for  intensity  stability.   Cavity  length 
tuning  allowed  vibration-rotation  line  selection,  with  a 
strong  hysteresis  effect.   A  computer  analysis  of  error 
signal  production  for  a  pyroelectric  detector  and  for  a 
HgCdTe  detector  indicated  that  phase  sensitive  detection 
of  the  first  harmonic  signal  component  gave  the  typical 
frequency  discriminant.   Operation  of  a  servo  loop  integrator 
for  the  feedback  loop  has  been  considered.   A  suitable  Stark 
Effect  absorption  cell  for  feedback  signal  generation  has 
been  constructed,  but  has  not  been  applied  to  the  laser. 
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I.   INTRODUCTION 

The  Sylvania  Model  941  CO   laser  was  purchased  to  do 
atmospheric  propagation  studies.   It  was  determined  then  to 
be  too  unstable  for  such  studies.   A  program  was  undertaken 
to  study  the  instabilities  and  to  stabilize  the  laser  both 
in  frequency  and  intensity.   The  novel  method  of  Stark  cell 
stabilization  was  chosen  as  the  best  route  to  take. 

Stark  cell  stabilization,  being  relatively  new,  has  not 
come  into  commercial  production  but  has  been  demonstrated  as 
feasible  in  the  laboratory  and  used  on  a  few  prototype 
systems.   As  a  result  there  is  very  little  information  in 
the  literature  as  to  the  exact  details  of  the  technique. 
Particularly,  it  was  required  to  know  how  to  produce  an 
error  signal  and  which  detector  would  be  most  suitable  for 
the  stabilized  system. 

This  thesis  describes  the  preliminary  accumulation  of 
knowledge  required  to  stabilize  the  Sylvania  laser  using 
the  Stark  cell  technique. 


II.   THEORETICAL 

A.   ROTATIONAL -VIBRATIONAL  SPECTRA 

The  quantum  states  of  a  molecule  arise  from  three  types 
of  motion  :   electronic,  vibrational,  and  rotational.   The 
electronic  charge  cloud  surrounding  the  molecule  may  be 
represented  mathematically  by  wave  functions.   This  wave 
function  representation  results  in  quantized  energy  states 
just  as  with  a  single  atom.   Motions  or  changes  in  the 
electronic  charge  distribution  may  result  in  a  changing 
electric  dipole  moment  in  which  case  radiation  is  either 
absorbed  or  released.   Electronic  transitions  correspond  to 
frequencies  primarily  in  the  visible  and  ultraviolet 
portions  of  the  spectrum. 

The  atomic  nuclei  in  the  molecule  themselves  are  not 
stationary,  but  vibrate  about  some  equilibrium  position.   For 
each  electronic  level  there  may  be  several  modes  of  vibration 
which  are  also  quantized.   These  vibrational  levels  result 
in  a  splitting  of  the  electronic  levels.   Transitions  be- 
tween vibrational  levels  correspond  to  frequencies  in  the 
near  and  middle  infrared  portions  of  the  spectrum.   Vibra- 
tional transitions  which  result  in  a  changing  electric  dipole 
moment  occur  via  absorption  or  release  of  photons. 

The  molecule  as  a  whole  may  also  rotate  about  various 
axes.  This  rotational  motion  is  also  quantized  and  gives 
rise  to  a  splitting  of  each  vibrational  level  into  rotational 


sublevels.   Transitions  between  rotational  levels  of  the 
same  vibrational  state  correspond  to  far  infrared  and 
microwave  frequencies.   As  in  the  other  two  types  of  motion 
only  those  transitions  which  result  in  a  changing  electric 
dipole  moment  will  occur  with  absorption  or  release  of 
radiation  at  these  frequencies.   The  above  described  energy 
level  scheme  is  represented  by  Fig.  1. 

Carbon  dioxide  is  a  linear,  symmetric,  triatomic  molecule 
It  can  vibrate  in  three  different  modes  as  pictured  in  Fig.  2 
These  are:   v, ,  the  symmetric  stretch  mode  in  which  the  two 
oxygen  atoms  move  in  opposite  directions  along  the  longi- 
tudinal line  of  symmetry  while  the  carbon  atom  remains 
stationary;  v  ,  the  bending  mode  in  which  the  oxygen  atoms 
remain  stationary  and  the  carbon  atom  moves  along  a  line 
perpendicular  to  the  longitudinal  line  of  symmetry;  and  v^, 
the  asymmetric  stretch  mode  in  which  the  oxygen  atoms 
remain  stationary  and  the  carbon  atom  moves  along  the  longi- 
tudinal line  of  symmetry.   Only  the  v2  and  v^  modes  give 
rise  to  an  electric  dipole  moment  and  can  thus  couple 
directly  to  electromagnetic  radiation.   The  v,  motion 
appears  strongly  in  Raman  scattering. 

The  v?  mode  can  occur  equally  in  two  different 
orthogonal  directions.   Thus  the  v2  mode  has  a  twofold 
degeneracy.   The  quantum  number  I   is  a  measure  of  the 
angular  momentum  about  the  longitudinal  axis  as  a  result 
of  the  v  motion.   The  vibrational  quantum  states  are  then 
denoted  by  [v,  v^   v.,] . 
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The  vibrational  energy  associated  with  a  given  state  is 
approximated  to  first  order  by 

Ev  =  (vl  +  I3  hc^l  +  ^v2  +  I}  hc^2  +  (V3  +  I}  hc^3      (1) 

where  v1  =  1388.3  cm  ,  x>2   =   667.3  cm  ,  v-  s  2349.3  cm1 ,  the 
frequencies  in  wave  number  units  associated  with  the  v1 ,  v?, 
and  V-  modes,  respectively.   v1 ,  v   and  v^  in  the  above 
equation  are  the  vibrational  quantum  numbers  and  take  on 
values  0 ,  1 ,  2  ,  ...  . 

The  C02  energy  levels  of  principal  interest  in  laser 
action  and  their  corresponding  frequencies  are:   [00°1]  at 
2349.3  cm1,  [10°0,02°0];[  at  1388.3  cm1,  [10°0  ,  02°0  ]  n  at 
1285.5  cm1,  and  [O^O]  at  667.3  cm1.   The  [10°0,02°0]I  and 
[10°0,  02°0]II  representations  follow  the  notation  used  by 
Tyte  [Ref.  1]  and  signify  the  fact  that  the  [10°0]  and 
[02°0]  levels  are  in  Fermi  resonance.   The  main  transitions 
of  interest  and  these  four  levels  are  represented  in  Fig.  3. 
In  the  figure  the  £  and  n  notations  refer  to  the  electronic- 
vibrational  angular  momentum.   This  vibronic  angular 
momentum  K  takes  on  values  0,  1,  2  ...  corresponding  to 
X,  n,  A  ...  ,  respectively  and  is  the  vector  sum  of  the 
vibrational  angular  momentum  about  the  internuclear  axis 
and  the  electronic  orbital  angular  momentum  about  the 
internuclear  axis.   The  I,  n,  and  A  states  of  a  molecule  are 
analogous  to  the  s,  p,  and  d  states  of  an  atom.   The  g 
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subscript  denotes  that  the  vibronic  eigenfunction  remains 
unchanged  upon  reflection  at  the  center  of  symmetry,  i.e., 
the  eigenfunction  is  of  even  symmetry.   The  u  denotes  an 
eigenfunction  of  odd  symmetry.   The  plus  superscript  denotes 
that  the  eigenfunction  remains  unchanged  when  reflected  at 
any  plane  containing  all  of  the  nuclei  [Ref.  2]. 

The  C0?  molecule  can  rotate  about  any  axis  perpendicular 
to  the  longitudinal  axis.   The  energy  of  rotation  is  given  to 
first  order  by 


t2 
Ej   =  J(J  +  1)  ^    J  =  0,  1,  2,  ...  (2) 


where  I  is  the  moment  of  inertia  and  J  is  the  rotational 
quantum  number.   From  Eq.(2)it  is  observed  that  the  energy 
difference  between  two  adjacent  rotational  states  is 


AEJ  =  EJ+1  -  EJ  =  I"  (J  +  «  (3) 


where  J  is  the  quantum  number  of  the  lower  state.   Substituting 
into  (3)  the  equation 


*Ej  =  £   ,  W 


solving  for  1/X  the  reciprocal  wavelength,  and  taking  the 
difference  in  wave  numbers  between  the  spectral  lines  one 
arrives  at 
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A<X)  "  17  -  XT.  "  HHT  AJ  •  ^ 


From  (5)  it  is  clearly  observed  that  for  two  adjacent  lines, 

i.e.,  AJ  =  ±1,  the  frequency  difference  is  a  constant  if  I  is 

a  constant.   In  general,  however,  I  is  not  a  constant,  but 

changes  slightly  with  different  vibrational  states.   This 

is  the  reason  for  the  observed  slightly  nonuniform  spacing 

of  the  rotational  lines.   Herzberg  [Ref.  2]  gives  various 

values  of  the  rotational  constant  for  the  corresponding 

-fc  i 

vibrational  levels.   A  typical  value  for  -. — r—  is  0.3899  cm 

for  the  [01  0]  vibrational  level.   This  nonuniform  spacing 

also  shows  up  in  the  departure  from  parabolic  of  the  potential 

curve  in  Fig.  1. 

The  selection  rules  originate  from  a  calculation  of 

the  electric  dipole  moment  matrix  element  taken  between  the 

initial  and  final  states; 

pfi  =  |/>*  er^dxl   .  (6) 

Transitions  between  states  which  yield  the  above  integral 
zero  are  forbidden.   For  a  harmonic  oscillator  the  integral 
yields  the  selection  rule  Av  =  ±1.   However,  as  the  vibra- 
tional energy  increases  and/or  the  molecular  structure 
becomes  more  complex,  the  vibration  departs  from  simple 
harmonic  and  Av  =  ±1  cannot  be  expected  to  hold.   From  (6) 
also  comes  the  selection  rule  for  C02  that  a  u  state 
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must  go  to  a  g  state  and  vice  versa.   But  u  to  u  and  g  to  g 
transitions  are  forbidden.   This  is  clearly  seen  in  Fig.  3. 

In  general  the  total  energy  of  a  molecule  to  first 
order  approximation  is  the  sum  of  the  energies  due  to  the 
electronic,  vibrational,  and  rotational  states; 

E  =  Ee  +  Ev  +  EJ  <7> 


The  I.R.  transitions  of  CO-  are  of  the  vibration-rotation 
type.   This  is  represented  schematically  in  Fig.  4  for 
transitions  between  two  vibrational  levels,  v.  and  Vr.   The 
rotational  selection  rule  AJ  =  ±1  causes  the  spectrum  to  be 

divided  into  two  branches:   a  "P  branch"  where  AJ  =  +1  and  an 

+    + 
MR  branch"  where  AJ  =  -1.   For  CO   this  is  a  Z   -  Z     vibronic 

transition.   Thus  AJ  =  0  is  forbidden,  which  would  otherwise 

yield  a  Q  branch  at  the  band  center  frequency  v  .   The 

frequencies  for  the  P  branch  are  represented  by 


v  =  vn    -  -iJy-  J    J  =  2,    4,  6  ...  (8) 

°    4tt  I 

and  for  the  R  branch 


v  =  v   +  -X-  (J  +  1)    J  =  0,  2,  4  ...        (9) 
0    4tt  I 


E    -  E 
v.     v£ 

where  v   =  r the  band  center  frequency  and  J  is  the 
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rotational  quantum  number  in  the  lower  vibration  level. 
Equations  (8)  and  (9)  are  easily  derived  from  Eq.  (2). 

Note  the  values  which  the  quantum  number  J  may  take  on. 
C02  is  a  linear  symmetrical  (point  group  D^)  molecule. 
The  two  identical  oxygen  nuclei  in  the  C0~  molecule  have 
zero  spin  and  obey  Bose  statistics.   Thus  all  antisymmetric 
levels  have  zero  statistical  weight  and  the  odd  rotational 
levels  of  a  Z      state  are  missing  as  well  as  the  even 

+  +       + 

rotational  levels  of  a  £   state.   Consequently  the  Z   to  I 

u  ^     J  u     g 

transitions  for  the  10.4  urn  and  9.4  ym  bands  have  the  simple 
structure  illustrated  in  Fig.  4  consisting  of  an  R  branch 
and  a  P  branch  with  alternate  lines  missing.   A  particular 
vibration-rotation  line  is  denoted  by  P(J)  or  R(J)  where  J 
is  the  quantum  number  of  the  rotational  level  in  the  lower 
vibrational  state.   A  few  of  the  more  important  lines  in  the 
10.4  urn  band  along  with  their  wavelengths  are  given  in 
Table  1.   Of  these  lines  the  P(20)  transition  is  the  strongest 

B.   EXCITATION  MECHANISM 

The  excitation  of  C02  into  the  [00°1]  level  is  readily 
accomplished  through  transfer  of  vibrational  energy  from 
nitrogen  molecules  in  the  v  =  1  vibrational  level  of  their 
ground  electronic  state.   From  Fig.  3  it  is  seen  that  the 
v  =  1  excited  vibrational  level  of  N-  is  about  19  cm  below 
the  [00°1]  level  of  CO-.   The  kT  average  thermal  energy  of 
molecules  at  room  temperature  is  about  210  cm  .   Thus  these 
two  levels  are  in  almost  exact  coincidence  and  a  high  cross 
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section  exists  for  vibrational  energy  transfer  through  two- 
body  collisions.   This  may  be  represented  by  the  equation 

N*(v  =  1)  +  CO2(00°l)  -*■  N2(v  =  0)  +  CO2(00°l)  -  19  cm1     (10) 

* 

where  N-  represents  a  vibrationally  excited  N?  molecule  in 

the  ground  electronic  state. 

N2  is  easily  produced  by  a  discharge  in  low  pressure 
nitrogen.   The  homonuclear  molecule  N-  has  no  permanent 
electric  dipole  moment;  thus  av=ltov=0  transition  via 
electric  dipole  radiation  is  forbidden,  causing  N-(v  =  1) 
to  be  metas table.   The  downward  transition  for  ^(v  =  1) 
is  governed  almost  solely  by  collisions  with  other  molecules 
and  the  container  walls. 

It  should  be  noted  that  the  excitation  cross  section 
of  CO   to  the  [10°0,  02°0]  level  by  N2  is  very  small  as  the 
difference  in  energy  is  more  than  900  cm  .   This  fact  along 
with  the  shorter  lifetime  of  the  [10°0,  02°0]  level  allows  a 
population  inversion  in  CO-  to  be  achieved  easily.   Helium 
is  often  added  to  increase  the  relaxation  rate  of  the  lower 
level. 

C.   LINESHAPES 

Analysis  of  a  spectral  line  of  an  atomic  or  molecular 
transition  reveals  it  to  be  composed  not  of  a  single  frequency 
only,  but  that  the  line  has  a  finite  width  or  frequency  spread 
with  an  intensity  profile  which  is  either  Lorentzian  or 
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Gaussian  in  the  limit.   The  frequencies  form  a  continuum  from 
zero  to  infinity  with  the  maximum  intensity  at  v  ,  the 
fundamental  frequency  of  the  transition.   These  two  pro- 
files are  compared  in  Fig.  5. 

The  Lorentzian  or  natural  lineshape  is  characteristic 
of  homogeneous  broadening.   In  homogeneous  broadening,  the 
absorbing  or  emitting  molecules  are  indistinguishable  as 
they  all  absorb  or  emit  at  the  same  frequency.   The  three 
most  common  causes  of  homogeneous  broadening  are:   (1)  the 
finite  spontaneous  lifetime  of  an  excited  state,  (2)  colli- 
sions of  molecules  with  phonons ,  and  (3)  collisions  of 
molecules  with  molecules  in  a  high-pressure  gas  which 
effectively  decreases  the  lifetime. 

The  Gaussian  lineshape  is  characteristic  of  inhomo- 
geneous  broadening.   Inhomogeneous  broadening  occurs  when 
the  absorbing  or  emitting  molecules  are  distinguishable  as 
a  result  of  their  absorbing  or  emitting  at  different 
frequencies.   The  two  most  common  causes  of  inhomogeneous 
broadening  are:   (1)  impurities  and  variations  in  a  host 
crystal  affect  the  energy  levels  and  thus  change  the  transi- 
tion frequency  of  nearby  emitting  molecules  and  (2)  the  high 
velocities  encountered  in  hot  plasma  gases  cause  the 
frequencies  to  be  Doppler-shif ted.   Obviously  for  a  C02 
laser  pressure  broadening  will  predominate  at  high  gas 
pressures,  resulting  in  a  Lorentzian  lineshape,  and  Doppler. 
broadening  will  predominate  at  low  gas  pressures,  resulting 
in  a  Gaussian  lineshape. 
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Consider  cause  No.  1  of  homogeneous  broadening  above. 
It  is  well  known  from  Fourier  transform  theory  that  a  pure 
single  frequency  implies  a  transition  which  occurs  over  an 
infinite  amount  of  time  which,  of  course,  does  not  occur. 
The  Lorentzian  lineshape  may  be  derived  by  assuming  the 
emitter  is  a  damped  harmonic  oscillator  with  an  exponential 
decay  time  x  [Ref.  3]. 


e(t)  =  EQe 


t/x 


cos  2ttv  t 
o 


(11) 


The  Fourier  transform  of  (11)  is  equivalent  to  the  convo- 

-t/x 
lution  of  the  Fourier  transform  of  e     and  the  Fourier 

transform  of  E   cos  2ttv  t : 


F{e(t)}  =  F{e"t/T}  *  F{E   cos  2ttv  t} 


r/0Qe-t/T  cos  2TTvtdt  -  i/°°e"t/T  sin  2rrvtdt]  * 
o  o 


r  E  E 

_£  6(v  .  Vq)  +  _o  6(v  +  VqJ 


1/x 

-in 


i2iTv 
T 


~T~1 


(1/x)"  +  4tt  v     (1/t)   +  4tt"v 


_£  5(v  .  v  )  +   o  5(v  +  vq)J   (12) 


where  the  lower  limit  of  integration  is  taken  as  zero  rather 
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than  negative  infinity  to  correspond  to  the  start  of  the 

emission.   Since  negative  frequencies  are  being  ignored  in 

this  analysis,  the  delta  function  at  v  =  +v   in  (12)  may 

be  taken  with  height  E   and  convolved  with  the  transform  of 
-t/x 


F{e(t)}    = 


1/x 


i2TTV 


(1/t)2    +    4tt2v2         (1/t)2    +    4tt2v2_ 


*  [EQ  5(v  -  vQ)]   (13) 


Equation  (13)  yields: 


1/t)E  i2tt(v  -  v  )E 

F{e(t)}  =  2 ^ 7   "  2 Z"2— ^ 7     '    (14) 

(1/tK  +  4tt  (v  -  vjz    (1/t)Z  +  4tt  (v  -  v  r 


Taking  the  square  of  the  absolute  value  of  (14)  (which 
corresponds  to  the  power  spectral  density  in  communications 
theory)  yields : 


F{e(t)} 


?      ?  ?  2       2 

(1/t)^E        +    4tt    (v    -    vj    Ert 

*•     '       J         O  0  0 

2 2 2~~2 

[(1/t)Z    +    4ttZ(v    -    v   )V 


2 2 1 

(1/t)Z    +    4ttZ(v    -    vQr 


(15) 


Normalizing  (15)  according  to_^  g(v)dv  =  1  yields 
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gO)  =  2   2/T2 7        •  (16) 

(l/iT  +  4tt  (v  -  v  )Z 


This  is  the  Lorentzian  lineshape  in  terms  of  x.   Defining 
Av  as  the  full  width  at  half  maximum  the  Lorentzian  may  be 
found  in  terms  of  Av .   From  the  definition  it  follows  that 


g(v0)/2  =  g(vQ  +  |^-)  .  (17) 


Substituting  in  the  proper  values  of  (16)  into  (17)  and 
carrying  out  the  algebra  yields 


i  =  ttAv   .  (18) 


Equation  (18)  is  worthy  of  closer  examination  at  this 
point.   It  is  noted  that  Av  is  a  measure  of  the  frequency 
spread  or  finite  line  width.   x  is  the  transition  time.   As 
t  decreases  the  frequency  spread  Av  must  increase.   Only 
if  x  were  infinite  would  Av  be  zero  implying  a  pure  single 
frequency  transition.   This  is  consistent  with  the  remark 
made  earlier  concerning  transitions  occurring  in  an  infinite 
amount  of  time. 

Substitution  of  (18)  back  into  (16)  yields  the  Lorentzian 
lineshape  in  terms  of  Av ; 
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g(v) 


2ir 


Av 

rAv    x 2  f  » 

(—  )       +    (v    -    vQ) 


(19) 


From    (18)    and    (19)    the    following   is   observed   to   be    true 


Av    = 


TTT  TTg(VQ) 


(20) 


Equation  (20)  says  that  as  Av  increases  or  x  decreases  the 
intensity  at  v  must  decrease.   But  if  the  intensity 
decreases  at  v   it  also  decreases  at  all  other  frequencies 
within  a  region  near  v  ,  frequencies  at  which  a  laser  will 
normally  oscillate.   This  is  one  possible  explanation  for  the 
typically  observed  decrease  in  power  output  from  a  laser 
with  an  increase  in  temperature.   An  increase  in  temperature 
would  increase  the  pressure  without  increasing  the  number 
density  and  pressure  broaden  the  gain  profile  resulting  in  a 
decrease  in  output  at  the  cavity  oscillation  frequency. 

Reference  4  gives  the  time  between  collisions  in  a  gas 
at  pressure  p  for  a  molecule  of  mass  M  and  a  diameter  a  as 


pa 


1   >/MkT 
2  v  TW 


(21) 


The  universal  gas  law  is  given  as 
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p    =   NkT 


(22) 


where   N   is    the   number   density   of  molecules    in   a   gas. 

Substituting    (22)    into    (21)    and   defining   the    collision   cross 

2 
section   as   Q   =   ira      yields 


=    1      /mtT 
TP       W$  v  T5TT 


(23) 


Substituting  (23)  into  (20)  then  gives  the  full  width  at 
half  maximum  due  to  pressure  broadening 


Av 


7T 


16kT 


(24) 


Equation  (24)  clearly  shows  that  Av  increases  as  the  square 
root  of  the  temperature. 

Consider  now  cause  no.  2  of  inhomogeneous  broadening. 
The  Gaussian  or  Doppler  broadened  lineshape  is  easily 
derived  by  beginning  with  the  Maxwell  velocity  distribution 
function  in  one  dimension  which  is  the  probability  that  a 
molecule  of  mass  M  in  a  gas  at  average  temperature  T  will 
have  an  x  velocity  component  v 


♦(vxD  - 


•  M  11/2 

iwrj       exP 


■Mv 
x 


2  -, 


(25) 
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The  Doppler-shif ted  frequency  is  given  as 


v  =  v   +  —  v 
o    c   o 


(26) 


which  upon  rearranging  gives 


v  =  (v 


0 


(27) 


Since  v  is  also  a  random  variable  and  <j>(v  )  is  a  probability 

A. 

distribution  function,  a  change  of  variables  can  be  made 
according  to  the  following  formula  [Ref.  5]  from  probability 
theory 


gO)  =  <|>(vx)  |J(vx/v) 


(28) 


where  J(v  /v)  is  the  Jacobian  and  in  this  case  equals 

.A. 

dv  /dv.   Using  (25)  and  (27)  in  (28)  and  normalizing  one 

arrives  at  the  probability  that  a  molecular  oscillator  of 

mass  M  and  resonant  frequency  v  in  a  gas  at  temperature  T 
emits  a  photon  at  frequency  v 


g(v)  = 


r _Mj  1/2 

2ttkT 


exp 


-Mc    r      s2 
(v  -  vQ) 


L2kTv 


T 


(29) 


This  is  the  Gaussian  lineshape  function  for  inhomogeneous 
broadening  and  has  a  full  width  at  half  maximum  given  by 
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AvD=  2vq  /^|n  2  (30) 


Mc 


It  was  stated  at  the  beginning  of  this  section  that  the 
lineshape  is  either  Lorentzian  or  Gaussian  in  the  limit. 
In  reality  several  broadening  mechanisms  occur  simultaneous- 
ly including  Stark  broadening,  van  der  Waals  broadening, 
and  wall  collisions  to  name  a  few.   In  the  simplified  case 
where  two  different  broadening  mechanisms  yield  comparable 
linewidths ,  but  which  are  considerably  larger  than  all  other 
broadening  mechanisms,  and  the  two  mechanisms  are  statis- 
tically independent,  the  observed  lineshape  is  equal  to  the 
convolution  of  the  two  individual  lineshapes 

Sobs(v)  =  r§l(v  '  ^')g2(v')dv'  (31) 

If  one  of  the  lineshapes  is  Gaussian  due  to  Doppler  broaden- 
ing and  the  other  is  Lorentzian  due  to  pressure  broadening, 
Eq.  (31)  yields  a  Voigt  profile  [Ref.  6].   The  total 
linewidth  in  this  case  may  be  estimated  according  to  [Ref.  7] 


Av  *  [(Av  )2  +  (Avp)2]1/2  (32) 


Equation  (31)  was  numerically  calculated  on  a  Hewlett- 
Packard  9810A  calculator  for  the  Lorentzian  and  Gaussian 
profiles  in  Fig.  5.   The  normalized  result  is  also  plotted 
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in  Fig.  5.   All  three  profiles  were  plotted  by  the 
calculator  on  the  same  scale  for  comparison  purposes. 
The  program  for  the  Voigt  profile  is  given  in  the  appendix 
for  the  inquisitive  reader. 

D.   LONGITUDINAL  AND  TRANSVERSE  MODES 

The  natural  resonant  frequency  of  a  laser  oscillator 
cavity  is  given  by 


v 


q£m 


cos 


/7T 


(I   +  m  +  1) 


R. 


0(1 


L_ 
R. 


2nL 


(33) 


where  q  is  the  longitudinal  mode  number,  I   and  m  are 
transverse  mode  numbers,  R,  and  R?   are  the  radii  of  curva- 
ture of  the  cavity  mirrors,  c  is  the  speed  of  light,  n  the 
refractive  index,  and  L  is  the  mirror  separation.   It  is 
easily  verified  from  this  equation  that  the  frequency 
separation  between  two  adjacent  longitudinal  cavity  modes  is 
given  by 


q+i   q 


2nL 


(34) 


A  laser  which  has  gain  over  this  frequency  difference,  i.e., 
a  linewidth  greater  than  or  equal  to  the  value  of  (34)  ,  will 
support  more  than  one  longitudinal  mode  at  a  time.   Thus 
a  single  mode  device  should  have  a  short  cavity  length. 
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The  frequency  or  mode  separation  for  adjacent  transverse 
modes  is  given  by 


vJ2,+m+l  "  v£+m 


cos 


Ci  -  JHU  - 


L  .  n 
RTj 


2nL 


(35) 


This  quantity  is  considerably  less  than  that  of  (34) ,  and 
often  less  than  the  line  width.   The  problem  of  higher  order 
modes  oscillating  simultaneously  is  sometimes  remedied  by 
the  use  of  an  aperture.   Bloom  [Ref.  8]  states  that  the 
lowest  order  mode  will  predominate  for  laser  bore  diameters 
less  than  3.5xco  where  u>  is  the  spot  size  on  the  output 
mirror. 

The  spot  size  for  a  resonator  with  mirrors  of  unequal 
radii  of  curvature  is  given  by  [Ref.  9] 


CO. 


AR  2 


2    u   L 

R,  -  L  R,  +  R. 


(36) 


Dividing  the  numerator  and  denominator  on  the  right-hand  side 
of  (36)  by  R12  and  taking  the  limit  as  ^  ^  »  gives  the 
spot  size  when  the  output  mirror  is  an  optical  flat; 


-i4  -  m2<*2  -  ^ 


(37) 
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E.   FREQUENCY  PULLING 

The  passive  cavity  linewidth  for  R  equal  to  the  square 
root  of  the  product  of  the  two  mirror  reflectances  and  a 
equal  to  the  average  distributed  loss  constant  is  given  by 
[Ref.  3] 

c[a  -  *  fcn  R]       _ 
Av   =  - ~  c(-1    K)  (38) 

where  the  approximation  holds  for  a  =  0  and  R  *  1 .   The 
cavity  linewidth  is  generally  less  than  the  Doppler  broadened 
linewidth  and  results  in  hole  burning  in  the  gain  curve  of 
an  inhomogeneously  broadened  system  at  the  frequency  or 
frequencies  of  oscillation. 

In  general  the  passive  cavity  resonance  v   does  not 
coincide  with  the  molecular  line  center  frequency  v    .   As  a 
result  the  oscillation  will  occur  at  a  frequency  near  v 
but  slightly  shifted  toward  v  .   This  is  called  frequency 
pulling.   This  frequency  pulling  is  due  to  the  negative 
anomalous  dispersion  curve  which  results  in  a  shift  in  the 
index  of  refraction  near  the  molecular  resonance.   The 
frequency  of  oscillation  of  the  laser  for  a  high  Q  cavity 
is  given  by 

,Avc  (39) 

v  "   vq  "  (vq  "  VlT 

where  Av  is  the  total  molecular  linewidth  due  to  Doppler  and 
pressure  broadening. 
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Since  the  longitudinal  mode  number  q  is  immensely 
larger  than  the  transverse  mode  numbers  I   and  m,  Eq.  (33) 
may  be  simplified  to 


qc 


2nL 


(40) 


Substitution  of  (38)  and  (40)  into  (39)  gives  for  the 
frequency  of  oscillation 


v  =  V 


2nL 


qc  (1  -  R) 
4irn  L 


voc(l  -  R) 

2-rrnL 


1_ 

Av 


(41) 


Equation  (41)  is  an  important  result  as  it  reveals  the  fact 
that  the  laser  oscillation  frequency  may  be  tuned  simply  by 
changing  the  cavity  length  or  broadening  the  linewidth. 
This  is  of  fundamental  importance  in  stabilization  tech- 
niques utilizing  a  piezoelectric  driven  mirror  in  the  cavity 

F.   INSTABILITY  AND  CAUSES 

Taking  the  derivative  of  (41)  with  respect  to  length 
yields 


dv  _  -qc 

cTL~   0  t  I 
2nL 


2ri    D^    v„c(l  "  R) 
qc  (1  -  R)  _   °     „ 

L  2im  L         2imL 


1_ 

Av 


(42) 


Substituting  Av   and  v   from  (38)  and  (40)  back  into  (42) 

gives  the  change  in  frequency  dv  for  a  change  in  length  dL 

as 
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£>  =    q  +  l(2v   -vj^  (43) 


The  quantity  in  the  parentheses  of  (43)  is  approximately 
v  ;  thus  (43)  is  simplified  to 


j     -v    v   Av 

HI    L    L   Av  L44j 


Equation  (44)  shows  that  when  Av   <<  Av  the  change  in 
frequency  with  a  change  in  length  is  large.   As  the  cavity 
linewidth  increases  or  the  molecular  linewidth  decreases 
the  laser  oscillation  frequency  fluctuates  less  with 
corresponding  length  fluctuations.   Normally  the  cavity 
linewidth  is  much  less  than  the  molecular  linewidth; 
consequently  the  primary  sources  of  instability  are  those 
which  affect  the  optical  length  of  the  cavity. 

Equation  (44)  also  yields  the  following  approximate 
expression  [Ref.  10] 


dv  =  -dL  (45) 

v     L 


Since  ■* —  is  proportional  to  the  thermal  coefficient  of 
linear  expansion  a  times  the  change  in  temperature,  it 
follows  from  (45)  that  the  frequency  of  oscillation  fluctuates 
with  temperature  according  to 
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p-   =  -adT  (46) 


Equation  (46)  shows  that  an  increase  in  temperature  yields 

a  decrease  in  frequency.   Mocker  [Ref.  11]  gives  a  coefficient 

o 

of  linear  expansion  for  low  expansion  quartz  of  2  x  10   /°C. 
With  this  value  of  a  the  laser  coolant  must  be  held  to 

within  1°C  in  order  to  achieve  a  frequency  stability  of  one 

7 
part  in  5  x  10  . 

Bridges  and  Patel  [Ref.  12]  point  out  that  there  is  a 
narrowing  of  the  Doppler  linewidth  with  decrease  in 
temperature,  which  is  evident  from  (30),  and  consequently 
one  would  expect  better  frequency  discrimination  between 
competing  modes.   Thus  if  mode  competition  is  a  problem 
where  an  interf erometric  device  or  grating  cannot  be  in- 
serted within  the  cavity,  lowering  the  coolant  temperature 

may  be  the  solution.   It  is  also  stated  in  Ref.  12  that  the 

-3/2 
optical  gain  is  approximately  proportional  to  T  '  .   Thus 

lower  temperatures  yield  higher  output  powers  as  well  as 

increased  stability. 

Pressure  fluctuations  are  another  source  of  instability, 

though  insignificant  in  sealed  lasers.   According  to  Mocker 

[Ref.  13]  frequency  shifts  due  to  gas  pressure  changes  of 

5-8  MHz  per  torr  have  been  observed.   This  shift  is  toward 

the  red  with  increasing  pressure.   The  shift  is  due  to  a 

change  in  the  refractive  index  of  the  gas  which  effectively 
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changes  the  optical  length  of  the  cavity.   For  a  sealed 
cavity,  pressure  fluctuations  would  result  from  changes  in 
the  gas  temperature  especially  during  warm  up. 

Acoustic  vibration  is  a  major  source  of  instability. 
It  also  results  in  cavity  length  fluctuations.   It  is 
minimized  by  utilizing  internal  cavity  mirrors  and  isolating 
the  cavity  as  a  whole  from  the  support  members  with  acoustic 
absorbers  such  as  foam  rubber  pads.   The  whole  laser  may 
then  be  enclosed  in  an  anechoic  box.   Significant  vibration 
results  from  the  circulating  coolant  water  and  fluctuations 
in  the  tap  pressure.   This  is  minimized  by  maintaining  the 
flow  rate  as  low  as  possible.   Air  cooled  systems  and  motor 
driven  choppers  can  generate  copious  amounts  of  acoustic 
noise. 

Fluctuations  in  the  excitation  current  result  in  both 
frequency  and  amplitude  instability.   Mocker  [Ref.  13] 
observed  a  shift  of  500-900  kHz/mA  shift  toward  the  blue 
with  increasing  excitation.   The  amount  of  shift  decreases 
with  increasing  excitation.   The  amplitude  or  output  intensity 
is  observed  to  have  an  optimum  value  of  excitation  and 
decreases  with  a  change  from  this  value  in  either  direction 
[Ref.  12] .   A  laser  which  is  to  be  amplitude  and  frequency 
stabilized  must  have  a  well  regulated  and  filtered  power 
supply. 

The  current  dependence  of  the  oscillation  frequency  is 
probably  due  to  a  change  in  the  refractive  index  of  the  active 
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medium  as  a  result  of  density  changes  of  the  electrons, 
ions,  and  neutral  molecules.   The  index  of  refraction  as  a 
function  of  the  number  densities  of  electrons  (N  ),  ions  (N-), 
and  molecules  (N  )  is  given  by 


,    N  e2 

n      -    1    -   4  — y  +  2tt    I   a,N.    +    2iTa      N         .                          (47) 

o                  2   m_   #l2  v      i    i               mm                                ^      J 

me   w  i 


A  reduction  in  the  refractive  index  may  arise  from  the  more 
rapid  creation  of  components  such  as  CO  and  NO  which  have 
a  smaller  electric  permittivity  than  C02  and  thus  contribute 
to  an  increase  in  the  frequency  of  oscillation. 

The  total  fractional  change  in  frequency  is  equal  to 
the  sum  of  the  fractional  fluctuations  [Ref.  10], 


*1   =  ^  *   *k      ,  (48) 

v    n    L   *  y     J 


where  An  is  found  from  (47)  and  AL  from  (45)  and  (46) . 

G.   STABILIZATION  TECHNIQUES 

As  was  noted  earlier  the  output  amplitude  of  a  laser  is 
primarily  a  function  of  the  discharge  current.   Therefore, 
current  regulation  should  eliminate  a  major  portion  of  the 
amplitude  fluctuations.   It  is  also  noted  that  due  to  a 
frequency-dependent  gain  curve  the  output  amplitude  is 
somewhat  dependent  on  the  frequency  of  oscillation. 
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Consequently  frequency  stabilization  will  also  produce 
some  amplitude  stability.  The  following  discussion  of 
some  of  the  techniques  used  focuses  on  frequency  stabilization 

In  the  above  discussion  of  instabilities,  it  was 
observed  that  each  fluctuation  effected  the  optical  length 
of  the  cavity,  thereby  changing  the  frequency.   It  follows 
then  that  frequency  stability  would  be  achieved  by  cavity 
length  correction.   This  is  often  accomplished  by  mounting 
one  or  both  morrors  of  the  cavity  on  a  piezoelectric  crystal 
which  expands  or  contracts  in  proportion  to  a  voltage 
applied  across  it.   An  error  signal  is  then  produced  by 
some  method  and  drives  the  piezoelectric  crystal.   A 
suitable  error  signal  must  meet  two  basic  criteria;  it  must 
be  of  a  magnitude  proportional  to  the  amount  of  length 
correction  required  and  it  must  indicate  the  direction  of 
the  required  length  correction. 

The  various  stabilization  techniques  may  be  basically 
divided  into  two  categories,  those  which  stabilize  to  a 
passive  cavity  resonance  and  those  which  stabilize  to  an 
atomic  or  molecular  transition  [Ref .  14] . 

An  example  of  stabilizing  to  a  passive  cavity  is  to  pass 
the  laser  beam  through  a  Fabry-Perot  etalon,  the  plate 
separation  of  which  is  modulated  at  an  audio  frequency. 
The  effect  of  the  audio  frequency  is  to  dither  the  cavity 
resonance  frequency  about  the  desired  optical  frequency, 
thus  changing  the  transmission  of  the  etalon.   The  output 
beam  of  the  etalon  is  a  phase  and  amplitude  modulated  signal 
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which  is  then  converted  to  an  electrical  signal  with  a 
photodetector.   The  signal  from  the  detector  is  compared 
with  the  audio  modulating  frequency  in  a  lock-in  amplifier. 
The  lock-in  amplifier  D.C.  output  is  the  error  signal. 
The  disadvantage  of  using  a  passive  cavity  is  that  it  is 
still  subject  to  thermal  and  acoustic  fluctuations. 

One  example  of  stabilizing  to  a  molecular  resonance 
is  to  use  the  Lamb  dip  of  the  Doppler  broadened  lineshape. 
The  Lamb  dip  or  power  dip  is  the  result  of  hole  burning 
in  the  gain  curve.   The  laser  frequency  may  be  dithered 
about  the  dip  and  the  laser  output  compared  with  the 
dithering  signal  to  produce  an  error  signal.   This  method 
has  the  inherent  disadvantage  of  having  a  very  narrow 
frequency  region  in  which  stabilization  can  occur.   If  the 
laser  frequency  should  inadvertently  drift  out  of  the  region 
of  the  power  dip,  the  error  signal  would  have  the  wrong 
polarity  and  force  the  laser  frequency  further  from  the 
desired  resonance. 

A  better  method  of  stabilizing  to  a  molecular  resonance 
is  to  use  an  external  gas-filled  cell  whose  center  absorption 
frequency  is  dithered  about  a  point  near  the  laser  frequency 
by  application  of  an  oscillating  magnetic  or  electric  field. 
This  type  of  center  frequency  modulation  utilizes  the  Zeeman 
effect  or  Stark  effect,  respectively.   The  output  of  the  cell 
is  then  phase  sensitively  detected  as  before  to  produce  an 
error  signal. 
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H.   STARK  CELL  STABILIZATION 

If  an  absorbing  or  radiating  molecule  is  immersed  in 
a  strong  electric  field  its  energy  levels  and  consequently 
its  spectrum  will  be  changed.   This  is  the  Stark  effect  and 
is  similar  to  the  Zeeman  effect  where  a  magnetic  field  is 
used  as  the  perturbation  rather  than  an  electric  field. 
The  application  of  the  electric  field  removes  the  degeneracy 
of  the  vibrational-rotational  levels  and  results  in  a  splitting 
of  the  levels  into  several  other  levels.   Kelly,  Francke ,  and 
Feld  [Ref.  15]  give  the  relevant  energy  level  diagram  for 
absorption  of  several  CO-  laser  lines  by  NH-D  and  show  the 
splitting  of  the  vibrational-rotational  levels. 

There  are  many  gases  whose  absorption  peaks  can  be 
brought  into  coincidence  with  C0~  laser  lines  by  the  Stark 
effect.   Martin,  Corcoran,  and  Smith  [Ref.  16]  give  a  table 
of  different  gases  and  the  laser  lines  absorbed.   The 
remainder  of  the  discussion,  however,  will  focus  on  the 
absorption  of  the  P(20)  CO   laser  line  in  the  10.4  ym  band 
by  NH-D.   Brewer,  Kelly,  and  Javan  [Ref.  17]  found  that  this 
coincidence  could  be  obtained  at  field  strengths  of  3298 
V/cm,  4404  V/cm,  and  6604  V/cm. 

The  D.C.  Stark  effect  can  be  represented  from  perturba- 
tion theory  to  first  order  by 


E  ■  E°  -  <$    |p  •  e  U  >  (49) 

m    m    ^m '  *    s '  ym 
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where  E°  is  the  unperturbed  energy  of  state  m,  \\>      is  the 


m 


unperturbed  wave  function,  p  the  molecular  dipole  moment, 
and  es  the  Stark  field  operator  [Ref.  19]. 

The  transition  frequency  between  two  Stark  shifted 
levels  is 


E   -  E 
mn      H      *  ^uj 


Substituting  a  time -dependent  Stark  field  of  the  form 
£s(t)  =  es  cos  oat  into  (49)  yields 


EmM  "  Em  "  ^m'P  '  eSIV  C0S  wt  ^51^ 


where  the  cosine  has  been  brought  outside  the  integral 

since  the  integration  is  over  distance.   It  is  also  assumed 

in  (51)  that  to  <<  vm  ;  otherwise  the  energy  levels  could 

mn 

not  follow  the  field.   Equation  (50)  then  becomes 


EmM    -    En(t) 

V    (t)   =  — r 

mnv  J  h 


EQ    -    E°    +    (<^    ip»eJ^>    -    <ty    |p-e    U   >)cos   cot 
ran  rnir      s|rn  rmir      s ' rm   '  rc_. 


Equation    (52)    can  be   rewritten   as 


v      (t)    =   v      +   Av '    cos   ojt    .  (53) 

mn  o 
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Equation  (53)  is  an  important  result  which  w^.11  be  used  in 
a  later  analysis.   It  shows  that  the  application  of  an 
oscillating  field  results  in  a  frequency  modulation  about 
a  center  frequency  v  .   Claspy  and  Pao  [Ref.  18]  solve  the 
equation  of  motion  using  the  density  matrix  approach  and 
show  that  for  Stark  field  frequencies  less  than  the  homo- 
geneous linewidth  Av,  the  absorption  coefficient  is  a 
single  Lorentzian  line  with  a  time  varying  center  frequency 
moving  synchronously  with  the  Stark  field.   The  result  is 
similar  to  (53) . 

Nussmeier  and  Abrams  [Ref.  19]  give  the  frequency  dif- 
ference between  the  P(20)  10.6  ym  line  center  and  the 
absorption  center  of  NH~D  as  a  function  of  applied  field  as 

5  =  2042  +  0.143 |M|E  (54) 

where  5  is  in  MHz,  M  is  the  z  component  of  the  rotational 
quantum  number  J  and  can  take  on  integer  values  from  -4  to 
4,  and  E  is  the  applied  electric  field  in  V/cm.   Equation 
C54)  shows  that  for  M  =  ±4  an  electric  field  of  3570  V/cm 
is  required  to  bring  the  P(20)  line  into  resonance  with 
NH2D.   Equation  (54)  differs  only  slightly  from  the  results 
found  by  Brewer,  Kelly,  and  Javan  [Ref.  17].   The  quantity 
0.143 | M |  in  (54)  is  called  the  Stark  tuning  rate  and  is 
reported  only  slightly  differently  in  the  literature.   For 
example,  Tangonan  and  Abrams  [Ref.  20]  report  a  tuning 
rate  of  NH2D  and  the  P(20)  line  of  .564  MHz  cm/V  for  M  =  4 
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which  gives  0.141  |M|,  Johnston  and  Melville  [Ref.  21] 
report  a  tuning  rate  of  .568  MHz  cm/v. 

A  proposed  stabilization  scheme  is  shown  in  Fig.  6, 
which  is  a  slightly  modified  version  of  that  given  in 
Ref.  19.   The  Stark  cell  D.C.  power  supply  biases  one  plate 
of  the  cell  to  bring  the  P(20)  line  into  resonance  with 
the  NH2D  gas.   This  occurs  at  a  field  strength  of  3570  V/cm. 
The  other  plate  is  driven  with  an  AC  signal  producing  a 
field  strength  of  approximately  150  V/cm  at  about  10  KHz 
[Ref.  19].   The  cell  is  filled  with  a  1:1  mixture  of  ND..  and 
NH3  at  approximately  2  torr.   The  exchange  occurs  rapidly 
to  produce  NH2D  [Ref.  21]. 

A  beam  splitter  reduces  the  radiant  flux  incident  on 
the  Stark  cell  below  lW/cm   to  avoid  saturation  effects. 
The  transmitted  intensity  is  then  detected  with  an  SBN  pyro- 
electric  detector  or  HgCdTe  detector.   The  detector  output 
is  fed  to  the  "signal  in"  on  the  lock-in  amplifier.   The 
audio  oscillator  signal  is  fed  to  the  "reference  in"  on 
the  lock-in  amplifier.   The  lock-in  determines  the  phase 
of  the  detector  signal  with  respect  to  the  audio  oscillator 
input.   The  polarity  of  the  D.C.  lock-in  output  signal  is 
determined  by  the  phase  of  the  detector  signal.   The  ampli- 
tude of  the  D.C.  output  is  proportional  to  the  amplitude 
of  the  detector  signal.   The  D.C.  lock-in  output  signal 
then  is  the  error  signal. 

The  error  signal  is  integrated  over  time  in  the  integrator 
The  integrator  output  drives  the  programmable  D.C.  power 
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supply  which  in  turn  drives  the  piezoelectric-mounted  laser 
cavity  mirror.   Thus  the  loop  is  completed.   The  exact 
details  of  error  signal  production  and  the  integrator  are 
covered  in  the  experimental  part  as  a  result  of  a  computer 
analysis . 
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III.   EXPERIMENTAL 

The  following  discussion  of  experimental  procedure  and 
results  is  specific  to  the  GTE  Sylvania  Model  941  CC>  Laser 
whose  characteristics  are  given  in  Appendix  A  and  to  the 
Stark  cell  whose  characteristics  are  given  in  Appendix  B. 
Unless  otherwise  indicated,  wherever  the  phrases  "the  laser" 
or  "the  Stark  cell"  appear,  it  is  understood  that  they 
refer  to  the  laser  and  Stark  cell  described  in  Appendices 
A  and  B. 

A.   LINESHAPES 

An  attempt  to  determine  the  lineshapes  of  the  laser 
transitions  experimentally  was  made  by  plotting  relative 
intensity  and  mode  as  a  function  of  piezoelectric  voltage. 
The  experimental  arrangement  is  schematically  represented 
in  Fig.  7.   The  x-axis  of  the  x-y  plotter  sampled  the 
piezoelectric  voltage  from  terminals  8  and  10  of  the  PZ 
power  supply  through  a  voltage  divider.   The  PZ  power  supply 
is  a  Kepco  Model  ABC  regulated  D.C.  supply  with  a  range  of 
0  to  1500  volts.   It  was  operated  in  the  voltage  program- 
mable mode  with  a  positive  ground.   The  programming  voltage 
was  a  regulated  0  to  25  volt  D.C.  supply.   The  battery  box 
and  reversing  switch  were  added  to  obtain  the  proper  biasing 
for  tuning  over  the  full  range  of  the  Kepco  with  a  minimum 
of  switching. 
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The  y-axis  of  the  x-y  plotter  sampled  the  relative 
intensity  of  the  chopped  laser  output  via  a  pyroelectric 
detector,  preamplifier,  and  a  filtered  rectifier.   The 
longitudinal  mode  or  vibrational-rotational  transition  was 
continuously  monitored  using  an  Optical  Engineering  C02 
Spectrum  Analyzer.   Frequency  could  not  be  monitored  due 
to  the  lack  of  precision,  but  the  longitudinal  mode  could 
be  identified.   Typical  results  are  given  in  Fig.  8  for  an 
increasing  PZ  voltage  and  in  Fig.  9  for  a  decreasing  PZ 
voltage. 

It  is  noted  that  the  plots  are  not  identical  but  a 
hysteresis  effect  exists  [Ref.  22].   The  P(20)  line  is 
shifted  in  the  direction  of  the  voltage  scan.   Various  scan 
rates  were  tried  in  order  to  ascertain  that  the  effect  was 
not  due  to  the  response  of  the  instruments.   Nine  to  ten 
different  lines  were  observed  over  the  tuning  range.   It 
is  also  noted  from  Figs.  8  and  9  that  the  series  of  lines 
begins  to  repeat  itself  after  approximately  1.4  kV. .  Since 
the  resonance  condition  requires  the  cavity  length  be  an 
integer  number  of  half  wavelengths,  this  voltage  change 


corresponds  to  a  length  change  of  half  a  wavelength  and  the 

AL 
Av 


piezoelectric  -r—  parameter  may  be  estimated  as 


AL  3  1/2(10.6  um)  .       f 

Av      1 .  4  kV 


46 


O 


O 

o 


3 
U. 

CO  UJ 
<  O 

< 
LU  f— 

Q_J 

OO 

^  > 

Qrsj 


> 


o 


CO 


CO 

< 

"Z.  LU 

lu  cr 
i-  o 


T 
< 

o 
en 


N 


> 


..  O 


CD 
LU 
D 


••  m 


Hb-O 


■■o 


ill 
o 
< 

_l 
o 
> 

N 

a 


A1ISN31NI    3AllV"13a 


47 


O 


o 


< 

o 

CD 
N 

a. 

t 


h- 

o 

if! 

k: 

CD 

_) 

< 

Ll 

h- 

_1 

CO 

o 

< 

> 

UJ 

a 
o 

IN 

CL 

^ 

O 

X 

ZT 

\- 

CO 

<r 

U) 

hi 

^L 

re. 

LU 

o 

h- 

LjJ 

^ 

Q 

> 
.10 


..O 


UJ 

o 


■»m 


"O 


UJ 

< 

-J 
O 

> 

N 


A1ISN31NI     3AI1V13H 


48 


It  is  noted  that  this  value  corresponds  with  the  value 
given  by  Sylvania.   This  calculation  assumes  that  the 
piezoelectric  crystal  has  a  linear  response  to  voltage. 

From  (55)  and  Figs.  8  and  9  the  apparent  linewidth 
of  the  P(20)  transition  is  estimated  to  be  55  MHz  which 
is  close  to  the  Doppler  linewidth  given  by  Sylvania.   The 
cavity  linewidth  calculated  from  (38)  with  a  =  0  is  19  MHz. 
The  pressure  broadened  linewidth  is  100  MHz  and  the 
Doppler  linewidth  is  50  MHz.   The  total  linewidth  using 
(32)  is  112  MHz.   Taking  the  differential  of  (59)  and  using 
19  MHz  for  Avc,  112  MHz  for  Av ,  and  55  MHz  for  the  change 
in  passive  cavity  resonance  dv  ,  gives  a  change  in  oscilla- 
tion  frequency  dv  of  46  MHz,  which  would  be  closer  to 
the  actual  linewidth  if  the  above  numbers  are  correct. 

It  is  concluded  that  the  PZ  voltage  may  be  used  to 
estimate  the  laser  oscillation  frequency;  however,  the 
actual  frequency  change  is  only  slightly  less  than  the  change 
indicated  by  the  PZ  voltage  due  to  the  frequency  pulling 
effect.   A  more  precise  means  of  measuring  the  actual 
frequency  of  oscillation  would  be  required  to  verify  (39) 
and  determine  the  average  distributed  loss  constant  a. 
For  gross  changes  in  the  frequency  of  oscillation,  using  the 
PZ  voltage  alone  will  result  in  insignificant  error.   To 
determine  small  changes  such  as  linewidths  a  more  precise 
method  of  measuring  frequency  is  required.   It  would  also  be 
necessary  to  insert  a  prism  or  grating  in  the  cavity  to 
prohibit  the  neighboring  modes  from  oscillating  until  the 
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half  power  point  was  reached.   From  Figs.  8  and  9  it  is 
observed  that  the  half  power  points  are  never  reached  on  the 
P(20)  line  before  the  neighboring  line  begins  to  oscillate. 
There  is  considerable  overlapping  of  the  lines.   Mode  compe- 
tition is  also  quite  prevalent  in  a  CO-  laser.   Thus  it  is 
impossible  to  measure  linewidths  in  the  manner  attempted  above 

B.   MODES  OF  OPERATION 

The  laser's  cavity  length  is  .45  m.   This  gives  from  (34) 
a  mode  separation  of  333  MHz  which  is  considerably  larger 
than  the  total  linewidth  of  112  MHz.   Consequently  the  laser 
should  normally  oscillate  in  a  single  longitudinal  mode, 
which  in  fact  it  did.   Occasionally  two  lines  were  observed 
to  oscillate  simultaneously  and  on  one  occasion  four  lines 
oscillated  simultaneously  at  a  coolant  temperature  of  about 
18°C;  the  P(22),  R(16)   R(24) ,  and  R(30)  lines.   Also  a  very 
faint  P(14)  line  of  the  9 . 4  urn  band  was  observed  to  oscillate 
simultaneously  with  P(22)  on  the  10.4ymband  at  19-20°C 
coolant  temperature  out. 

The  laser  utilizes  a  dielectric  flat  for  the  output 
mirror  and  a  3  m  totally  reflecting  mirror  mounted  on  a 
piezoelectric  crystal.   Using  (55)  gives  a  transverse  mode 
separation  of  42  MHz,  which  is  less  than  the  total  line- 
width.   From  (37)  the  spot  size  on  the  output  mirror  is 
1.9  mm.   The  laser's  bore  diameter  is  6  mm.   Thus  from  the 
discussion  preceding  Eq.  (36)  the  laser  is  only  marginally 
diffraction  limited  to  the  lowest  order  transverse  mode.   On 
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one  occasion  the  TEM10/TEMQ1,  the  "doughnut"  mode,  was 
observed.   On  another  occasion  efforts  to  reproduce  this 
mode  yielded  only  the  TEM   mode  throughout  the  full 
range  of  the  PZ  voltages. 

C.   INSTABILITIES 

Experimentation  began  with  minimization  of  60  Hz 
noise.   The  voltages  on  the  ground  terminals  of  all 
115  VAC  receptacles  in  the  lab  were  measured  with  respect 
to  a  drainpipe  ground.   One  receptacle  ground  was  found  to 
have  a  52  mV  120  Hz  signal.   Its  ground  strap  was  isolated 
from  the  grounds  of  the  other  receptacles  in  order  to 
minimize  the  signal  on  them.   The  receptacle  with  the 
lowest  ground  signal  (5  mV  peak  to  peak  at  180  Hz)  was 
used  as  the  source  for  all  measuring  instruments  and  power 
supplies.   Noisy  receptacles  were  used  for  vacuum  equipment, 
chopper,  and  heat  lamps.   Ground  loop  signals  were  eliminated 
only  by  trial  and  error. 

The  laser  power  supply,  the  GTE  Sylvania  model  941P, 
is  unregulated  and  minimally  filtered.   It  was  found  to  have 
a  one  to  ten  per  cent  current  ripple  depending  on  the  dis- 
charge current.   This  was  entirely  unsuitable  for  stabiliza- 
tion of  the  laser  as  the  current  ripple  showed  up  in  an 
amplitude  modulated  laser  output.   An  L-section  filter 
consisting  of  a  .5  uf  capacitor  and  a  45  H  choke  (which  was 
actually  the  secondary  of  a  13  kV  transformer)  was  assembled 
and  connected  to  the  high  voltage  output  at  the  laser  power 
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supply.   Details  of  the  filter  are  given  in  Appendix  E.   The 
filter  reduced  the  current  ripple  and  consequently  the  laser 
output  fluctuation  by  a  factor  of  50  to  approximately  .1%. 
Since  a  101  Stark  cell  modulation  was  expected,  this  filter 
was  considered  to  be  suitable  for  short-term  stability 
purposes.   For  long-term  stability,  a  well  regulated  power 
supply  was  still  required.   Plots  of  laser  intensity  as  a 
function  of  excitation  current  were  made  using  an  x-y 
plotter.   The  x-axis  input  was  the  voltage  drop  produced 
by  the  excitation  current  across  a  small  resistance  element. 
Typical  results  are  given  in  Fig.  10.   The  relative  intensity 
variation  is  approximately  14%. 

The  experimental  setup  for  determining  the  effect  of 
cooling  water  temperature  on  the  intensity  and  longitudinal 
mode  is  schematically  represented  in  Fig.  11.   The  coolant 
output  temperature  was  sampled  using  a  copper-constantan 
thermocouple  referenced  to  an  ice  bath.   The  voltage  drop 
across  the  thermocouple  was  measured  with  a  nanovoltmeter 
on  the  3  millivolt  scale.  The  output  of  the  nanovoltmeter 
drove  the  x-axis  of  an  x-y  plotter.   The  y-axis  sampled 
the  chopped  laser  intensity  via  a  pyroelectric  detector, 
preamplifier  and  filtered  rectifier.   The  coolant  tempera- 
ture was  slowly  increased  over  a  period  of  about  one  hour 
from  7°  to  40°C  by  circulating  the  coolant  through  an  ice 
bath  and  warming  the  ice  bath  with  a  heating  element.   The 
longitudinal  mode  was  continuously  monitored  using  a  spectrum 
analyzer.   Mode  changes  were  noted  and  recorded  on  the  plot 
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after  each  run.   A  table  of  temperature  as  a  function  of 
thermocouple  voltage  is  given  in  Appendix  F. 

Figure  12  is  a  plot  with  the  PZ  voltage  held  constant. 
Note  that  the  series  of  lines  from  P(18)  to  R(22)  is  repeated 
in  identical  order;  however,  there  is  little  similarity  to 
the  intensity  profiles.   The  P(18)  line  is  repeated  twice 
over  a  temperature  range  of  8.8  to  37.9°C.   Since  the 
resonance  condition  requires  the  cavity  length  to  be 
exactly  an  integer  number  of  half  wavelengths  of  the  oscilla- 
tion frequency,  this  temperature  range  corresponds  to  a 
thermal  expansion  of  two  half  wavelengths  or  one  wavelength 
of  the  P(18)  transition.   Thus  the  thermal  coefficient  may 
be  determined  from  Eqs .  (45)  and  (46)  as 

a  =  I  dL  =   1   (10-57  x  10'6)  _     in-7/r°         rw 
a   L  aT   TUT  29.1UC -  8  x  10   /C   .       (56) 

Note  that  this  calculation  assumes  that  a  temperature  change 
in  coolant  out  corresponds  to  an  equal  temperature  change 
in  the  laser  tube  wall.   See  Appendix  G  for  proof  of  this. 
The  plot  shows  an  overall  decrease  in  intensity  output 
with  increasing  temperature,  by  a  factor  of  about  two  in  a 
temperature  range  of  33C°.   It  is  also  obvious  that  the  P(20) 
transition  is  the  strongest  in  the  group.   P(22),  R(16)  and 
R(24)  form  a  competing  trio.   Mode  competition  is  almost 
always  present  between  these  lines.   All  three  or  pairs  of 
these  lines  are  often  observed  oscillating  simultaneously. 
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It  would  therefore  be  wise  to  avoid  stabilizing  on  any  one 
of  these  three  lines. 

Figure  13  shows  the  maximum  intensity  of  the  P(20)  line 
as  a  function  of  temperature.   The  intensity  decreases 
almost  linearly  with  increasing  temperature.   The  region 
between  7  to  27°c  shows  mode  competition.   It  is  also  noted 
from  Fig.  13  that  as  temperature  increases  or  the  cavity 
expands,  the  PZ  voltage  had  to  be  increased  in  order  to 
maintain  the  P(20)  line.   It  follows  that  an  increasing  PZ 
voltage  corresponds  to  a  decreasing  cavity  length  and  the 
sign  on  the  ^y  parameter  given  by  (55)  should  be  negative. 

An  experiment  was  performed  to  determine  the  effect 
of  coolant  flow  rate  fluctuations  on  the  laser  cavity. 
The  flow  rate  was  adjusted  within  the  limits  set  by  Sylvania 
from  11  ml/s  to  30  ml/s  at  a  constant  temperature  input. 
The  temperature  output  changed  as  expected  from  23.8  to 
22.8°C.   Using  (56)  this  one-degree  temperature  change 
corresponds  to  a  length  change  of  .36  um.   In  this  experiment 
the  PZ  voltage  was  adjusted  to  maintain  maximum  intensity 
on  the  P(20)  line,  and  was  adjusted  from  .42  kV  to  .24  kV. 
Using  (55)  this  voltage  change  corresponds  to  a  length  change 
of  .68  um  which  is  almost  twice  the  length  change  accountable 
for  temperature.   It  follows  that  a  flow  rate  fluctuation  of 
19  ml/s  results  in  a  length  fluctuation  of  approximately 
.32  um.   From  (45)  this  corresponds  to  a  frequency  fluctua- 
tion of  2  2.6  MHz. 
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It  is  concluded  from  these  experiments  that  temperature 
and  flow  rate  fluctuations  are  the  dominant  causes  of 
instability  in  the.  laser.   Current  fluctuations  can  be 
easily  remedied  with  a  filtered  current  regulated  power 
supply.   If  a  circulation  pump  is  used  for  cooling,  the 
system  should  be  well  designed  to  minimize  pressure  and 
flow  rate  fluctuations.   Experiments  could  not  be  carried 
out  to  determine  the  effects  of  acoustic  vibration  due  to  a 
lack  of  a  method  for  precisely  measuring  small,  rapid 
fluctuations  in  frequency. 

D.   ERROR  SIGNAL  ANALYSIS  FOR  A  PYROELECTRIC  DETECTOR 

Drawing  an  analogy  with  electrical  communication  theory 
it  is  observed  that  the  stabilization  technique  is  similar 
to  frequency  demodulation.   The  input  laser  beam  is 
frequency  modulated  with  noise.   It  is  desired  to  have  some 
sort  of  frequency  discriminator  which  will  sense  the  shift 
in  frequency  and  produce  an  output  voltage  linearly  pro- 
portional to  the  shift.   This  is  precisely  what  FM  demodu- 
lation accomplishes  [Ref.  25].   The  stabilization  servo 
loop  is  completed  when  the  demodulation  information  is  fed 
back  to  the  laser  to  make  corrections  to  minimize  the  noise 
modulation. 

The  exact  details  of  producing  an  error  signal  were 
not  clear  in  the  literature.   Questions  which  needed  to  be 
resolved  were:   is  the  error  signal  produced  by  slope  or 
peak  detection,  how  does  dithering  produce  an  error  signal, 


59 


is  the  composite  detector  signal  or  one  of  its  harmonics 
monitored,  and  which  harmonic  should  be  monitored? 

In  order  to  search  for  a  suitable  error  signal  we 
begin  by  assuming  the  transmitted  intensity  of  the  cell 
of  length  L  is  of  the  form 

I(v)  =  IQ(v)  exp  [-y(v)L]  (57) 

where  IQ(v)  is  the  input  intensity,  y(v)  is  the  absorption 
coefficient,  and  I(v)  is  the  output  intensity  as  a  function 
of  frequency.   It  is  assumed  that  I(v)  is  a  narrow  spectral 
laser  line  of  half  width  much  less  than  the  cell  absorption 
linewidth.   Thus  IQ(v)  may  be  approximated  by  a  Dirac  delta 
function,  <5(x).   If  the  frequency  of  the  laser  is  allowed  to 
vary  randomly  with  time  over  the  linewidth  of  the  absorption 
cell,  and  the  linewidth  of  the  absorption  cell  is  assumed  to 
be  much  less  than  the  laser's  molecular  linewidth  so  that  the 
laser  intensity  does  not  change  appreciably  over  the  cell 
linewidth,  then  (57)  becomes 

Iv(t)  =  /°Io6(x(t)  -  v)e"Y(v)L  dv  =  IQ  exp  [-Y(x(t))L]      (58) 

where  x(t)  is  a  random  variable  of  time  and  takes  on  values 
of  frequency  within  the  cell's  absorption  linewidth.   The 
subscript  v  on  I  denotes  that  I  is  approximated  by  a  delta 
function,  i.e.,  all  the  intensity  is  concentrated  at  a 
single  frequency. 
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Equation  (58)  shows  that  the  frequency  modulation  x(t) 
has  been  converted  to  an  intensity  modulation  I(t)  by  simply 
passing  the  laser  beam  through  a  selective  absorber  or 
optical  notch  filter.   It  is  possible  to  derive  an  error 
signal  from  I(t)  by  stabilizing  about  a  point  on  the  slope 
of  y(v).   This  is  called  slope  detection,  analogous  to  slope 
detection  in  radio  FM  demodulation.   Two  problems  arise  with 
slope  detection;  spurious  intensity  variations  in  the  laser 
beam  also  show  up  in  the  error  signal  and  the  range  of  the 
linear  slope  is  quite  small  [Ref .  23] .   To  partially  over- 
come these  problems,  it  is  desirable  to  obtain  an  error 
signal  which  is  similar  to  the  frequency  to  voltage  charac- 
teristic of  a  common  FM  balanced  discriminator,  where  the 
voltage  goes  to  zero  at  the  center  frequency,  increases 
positively  on  one  side  of  the  center  frequency,  and  increases 
negatively  on  the  other  side. 

Consider  the  effect  of  modulating  the  center  frequency 
of  the  absorption  cell  with  an  audio  signal  of  frequency  u> 
according  to  (53) ,  and  using  phase-sensitive  detection  in 
order  to  determine  zero  time.   This  is  schematically  repre- 
sented in  Fig.  14.   The  modulating  signal  will  dither  the 
absorption  profile  about  v   at  frequency  w.   The  amplitude 
of  the  dither  is  Av ' .   It  is  assumed  that  the  dither  is 
much  more  rapid  than  the  random  variation  of  the  laser 
frequency  v.   As  will  become  more  apparent  later  the  reason 
for  the  assumption  is  to  detect  x(t) .   Only  fluctuations 
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slower  than  to  can  be  detected.   The  end  result  of  the 
analysis  to  follow  is  that  x(t)  amplitude  modulates  a 
"carrier  frequency"  co.   It  would  be  desirable  to  make  u> 
extremely  large  in  order  to  detect  as  much  noise  as  possible. 
However,  as  Claspy  and  Pao  [Ref.  18]  show,  the  upper  limit 
for  oj  is  the  homogeneous  linewidth.   This  is  reasonable  if 
one  considers  that  from  (18)  the  linewidth  was  derived  from 
the  decay  time.   It  is  now  observed  that  the  decay  time  is 
also  the  absorption  time.   Thus  if  a  molecule  does  not 
remain  in  resonance  with  the  passing  radiation  for  the 
duration  of  the  absorption  time,  the  molecule  will  not 
absorb.   In  other  words,  the  Stark  cell  is  bandwidth 
limited  by  the  homogeneous  linewidth. 

The  detector  used  in  this  part  of  the  analysis  has  an 
output  voltage  proportional  to  the  time  rate  of  change  of 
the  incident  intensity.   An  example  of  such  a  detector  is  a 
pyroelectric .   This  time-differentiating  detector  will  be 
referred  to  as  "detector  A"  in  this  analysis. 

The  absorption  coefficient  of  the  cell  is  assumed  to 
have  a  Lorentzian  shape  from  (19) . 


YM  =  Sg(v)  =  ^ (59) 


2,[(^)2  +  (V  -  vQ)2] 


where  S  is  the  line  strength  and  is  a  constant.   The  absorp 
tion  coefficient  is  given  a  time  dependence  by  oscillating 
the  center  resonant  frequency  according  to  (53). 
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Y(v,t)  = 


2tt 


SAv 

l-j— )   +  (v  -  v   -  Av  '  cos  oat) 


(60) 


The  difference  in  frequency  between  the  laser  beam  and  the 

average  center  resonant  frequency  v   is  defined  as 

5  =  v  -  v   ,  where  the  symbol  used  is  not  to  be  confused 

with  the  Dirac  delta  used  earlier.   Equation  (60)  then 

becomes 


Y(v,t)  = 


SAv 


2tt  ["(^-)2  +  (5  -  Av'  cos  cot)2 


(61) 


Substituting  (61)  into  (58)  the  intensity  out  as  a  function 
of  time  is  given  by 


I(t)  =  IQ  exp 


SAvL 


2tt 


£±)1    +  (o  -  Av'  cos  cot)2 


(62) 


Detector  A  output  as  a  function  of  time  is  a  voltage 
proportional  to 


£(t)  =  d\W    =  -I0SL  dg(dl't]    exp  {-SL  g(v,t)> 


(63) 


As  was  assumed  earlier  the  oscillation  frequency  is  much 
greater  than  the  random  variation  of  the  laser  frequency  v 


64 


Thus  6  is  approximately  constant  and  the  derivative  with 
respect  to  time  of  (62)  may  be  taken  holding  it  as  such. 
Taking  the  derivative  of  (59),  expanding  all  squared  terms 
and  using  common  trigonometric  identities  yields 


dg(v,t)  -SAvAv'u>(6    sin  cut    -   ^  sin.   2u)t) 

dt  •  ttLA+B-cos    cot   +   C'cos    2ojt   +   D-cos    3cot    +   E-cos    4oot] 


where 


A   =    54    +    352(Av')2    +   |(Av')4    +    (^)  2    262    +    (^i)2(Av')2    +    (^) 4 


B    =    -4S3Av'    -    36(Av')3    -    46Av»(^-)2 


2 


C    =    362(Av')2    +   ^r(Av')4    +    (Av')2(^i)2  (64) 

D   =    -6(Av')3 


E    =   J-(Av')4 


Equation  (63)  can  be  simplified  if  y(v,t)L  <<  1.   The 
exponential  then  is  approximately  one  and  f(t)  is  pro- 
portional to  (64)  .   Note  that  S  could  be  determined 
experimentally  by  applying  a  D.C.  field  and  measuring  the 
transmitted  intensity  when  v  =  v   and  when  |v  -  v  |  >>  Av . 

Equation  (64)  was  plotted  on  a  HP  9810A  calculator 
for  various  values  of  5  and  Av ' .   The  full  width  at  half 
maximum  Av  and  the  period  of  the  fundamental  oscillation  T 
were  arbitrarily  set  to  one  in  the  program;  thus  oo  the 
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fundamental  frequency  is  2tt .   5  is  varied  between  -1.0  and 
1.0  and  Av'  between  0.1  and  1.0.   The  results  are  plotted 
in  Figs.  15,  16,  17,  and  18.   Figure  15  is  for  5=0  and 
Av '  =  .1  and  .7.   Note  the  frequency  is  twice  the  funda- 
mental.  This  would  be  expected  for  a  laser  line  oscillating 
symmetrically  about  the  peak  of  a  transfer  characteristic. 
From  (64)  it  is  also  noted  that  when  6  =  0  no  fundamental 
or  first  harmonic  exists.   Figure  15  also  reveals  the  non- 
linearity  produced  by  overdriving  with  a  large  Av ' .   Though 
the  output  signal  is  greater,  higher  order  harmonics  are 
obviously  present,  producing  a  distorted  sin  2u>t. 

Figure  16  is  for  <5  =  .  2  and  Av '  =  .1  and  .7.   It  is 
obvious  that  some  first  harmonic  is  present.   It  is  noted 
that  when  Av '  >  5  the  dither  carries  the  line  to  both  sides 
of  the  y(v)  peak,  thus  producing  a  large  amount  of  second 
harmonic  but  because  the  dither  is  unsymmetrical  about 
the  peak  of  y(v)  a  small  amount  of  first  harmonic  is 
generated.   When  5  >  Av f  second  harmonic  is  virtually 
nonexistent. 

Figure  17  is  for  5  =  .5  and  Av '  =  .1  and  .7.   This  is 
the  region  of  slope  detection.   A  considerable  amount  of  first 
harmonic  is  obviously  present.   Very  little  second  harmonic 
is  present  even  when  Av '  is  greater  than  6.   Figure  18  shows 
a  180°  phase  shift  in  going  from  a  positive  6  to  a  negative  <5 . 

It  is  concluded  from  these  four  plots  that  the  basic 
requirements  for  generating  an  error  signal  are  present;  the 
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amplitude  of  the  signal  indicates  the  magnitude  of  5  and  the 

phase  of  the  signal  indicates  the  sign  of  5.   Thus  the 

distance  of  the  laser  frequency  from  the  reference  v   is 

known  and  the  direction  from  v   is  known.   It  is  still 

o 

desired,  however,  to  obtain  the  typical  FM  frequency  dis- 
criminant and  to  determine  which  harmonic  to  monitor  or 
if  the  composite  signal  should  be  monitored. 

It  is  revealing  to  decompose  the  signal  of  (64)  into 
all  the  harmonics  which  make  up  the  signal  and  examine 
each  harmonic  singly.   Any  periodic  signal  of  fundamental 
period  T  which  is  bounded  and  real  may  be  expanded  into  a 
Fourier  series  according  to 


rr^^        ao    r  t  nZiTt   u   _.   n2irt^ 

f  CO  =  —  +  I    (a  cos  -J —  +  bn  sin  -^ — ) 

n=l 


where 


an  =  1  /  f(t)  cos  S2lt  dt  ,  n  =  0,  1,  2 


T 

b  =  I  f   f(t)  sin  ^^  dt  ,  n  =  1,  2,  ...  (65) 

n   T  ^   v        T 
o 


For  the  complex  signal  given  by  (64)  the  integration 
cannot  be  performed  analytically;  however,  it  can  be  per- 
formed numerically  on  a  computer.   It  is  known  that  for  an 
odd  function  of  t,  i.e.,  f(-t)  =  -£ (t) ,  which  is  what  the 
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plots  show  (64)  to  be,  all  the  a  's  will  integrate  to  zero 
over  one  period.   Thus  only  the  b  's  need  be  calculated. 
For  a  fundamental  period  of  T  =  1,  Eq.  (65)  reduces  to 
solving  the  problem 


f(t)  =  I   b   sin  (n2irt) 
n=l 


bn  =  2  f   f(t)  sin  (n2TTt)  dt   n  =  1,  2,  ...      (66) 
o 


for  f (t)  proportional  to  (64) . 

A  program  for  the  HP  9810  was  written  to  do  the  numerical 
integration  of  (66)  up  to  n  =  10  for  various  values  of  5 
and  Av '  and  print  the  results.   It  was  concluded  that  the 

Fourier  coefficients  b   for  harmonics  greater  than  n  =  3 

n  ° 

were  too  small  to  be  of  any  value  in  stabilization.   A 
program  was  then  written  to  calculate  the  Fourier  coefficient 

b   for  a  particular  n  preset  into  the  calculator  for  200 

n       r  r 

different  values  of  6  from  -1.0  to  1.0  and  plot  the  results. 

As  before  Av  and  T  were  arbitrarily  set  to  one.   Av '  and 

o 

6  have  the  same  units  as  Av  and  are  to  be  viewed  as  a 
fraction  of  the  linewidth  Av .   Av  may  be  scaled  up  to  any 
desired  linewidth;  <5  and  Av '  are  then  scaled  up  with  the  same 
scaling  factor.   This  program  is  given  in  Appendix  D. 

Plots  of  the  first  three  Fourier  b  coefficients  of 
(64)  as  a  function  of  5,  with  S  -  1,  are  shown  in  Figs.  19, 


72 


20,  21,  and  22.   Figures  19  and  20  are  the  typical  shape 
of  an  FM  discriminant.   These  plots  represent  the  relative 
magnitude  of  the  D.C.  error  signal  produced  by  a  lock-in 
amplifier  when  only  a  single  harmonic  is  phase  sensitively 
detected  with  the  oscillator  signal  to  the  cell  as  the 
reference.   Other  harmonics  could  be  filtered  out  using  a 
narrow  band  pass  filter.   Figures  19  and  20  show  that  the 
first  harmonic  should  be  monitored.   Note  the  linear  region 
near  6=0.   When  the  laser  frequency  is  at  v   the  error 
signal  is  zero.   As  the  laser  frequency  drifts  to  the  right 
of  vq  for  a  positive  6  the  error  signal  increases  linear 
and  negative.   A  drift  to  the  left  of  v  would  give  a  posi- 
tive error  signal.   This  is  called  peak  detection  as  the 
laser  frequency  is  stabilized  about  the  peak  of  the  absorp- 
tion profile. 

Figure  20  shows  that  a  large  dither  amplitude  produces 
a  suitable  error  signal.   Note  that  the  peaks  are  further 
away  from  center  than  for  a  small  dither  amplitude.   Figure 
21  shows  that  the  second  harmonic  is  unsuitable  for  pro- 
ducing an  error  signal.   However,  it  might  be  used  in 
conjunction  with  the  first  harmonic  for  making  a  gross 
correction  to  the  laser  cavity.   Note  that  if  the  laser 
frequency  drifts  beyond  the  peaks  of  the  first  harmonic 
signal,  the  second  harmonic  goes  negative  indicating  that 
the  range  of  stability  has  been  exceeded  and  a  gross 
correction  to  the  laser  cavity  in  the  direction  indicated 
by  the  sign  of  the  first  harmonic  is  required  in  order  to 
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bring  the  laser  frequency  back  into  the  linear  portion  of 
the  discriminant. 

Figure  22  shows  that  the  third  harmonic  may  be  used 
to  produce  an  error  signal  if  the  dither  amplitude  is 
sufficiently  large. 

It  is  concluded  that  an  error  signal  may  be  produced 
by  phase  sensitive  detection  of  the  fundamental  or  first 
harmonic  of  the  detector  signal.   The  lock-in  amplifier 
establishes  zero  time  from  the  reference  signal.   The 
phase  of  the  detector  harmonic  will  either  be  0  or  180°. 
The  polarity  of  the  error  signal  is  determined  by  the 
phase  of  the  harmonic  being  monitored.   The  amplitude  of 
the  error  signal  is  proportional  to  the  amplitude  of  the 
harmonic,  and  for  the  first  and  third  harmonics  the  error 
signal  is  linearly  proportional  to  the  frequency  difference 
5  in  the  region  near  5=0. 

E.   ERROR  SIGNAL  ANALYSIS  FOR  A  HgCdTe  DETECTOR 

Consider  now  using  a  nondiff erentiating  detector.   The 
detector  has  an  output  voltage  proportional  to  the  absolute 
incident  intensity.   Such  a  detector  would  be  a  HgCdTe 
detector.   This  type  of  detector  will  be  referred  to  as 
"detector  B."  With  detector  B  in  the  loop  the  detector 
signal  will  be  proportional  to  I(t)  as  given  in  (58). 
Equation  (58)  can  be  approximated  to  first  order  according 
t0  ex  =  i  +  x  for  x  small.   Thus  the  detector  output  voltage 
is  proportional  to  (61)  and  should  contain  a  large  D.C. 
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term  due  to  expanding  the  exponential.   A  plot  of  (61)  as 

a  function  of  time  showed  that  this  was  indeed  the  case. 

The  function  is  also  an  even  function,  i.e.,  f(-t)  =  f (t) ; 

thus  expansion  is  a  Fourier  cosine  series  and  only  the  a  's 

J  n 

in  (65)  need  be  calculated. 

Plotting  the  coefficients  of  the  harmonics  as  a  function 
of  5  as  before  should  show  whether  or  not  detector  B 
yields  a  suitable  error  signal.   Plots  for  the  first, 
second,  and  third  harmonics  are  given  in  Figs.  23,  24,  and 
25.   The  results  are  similar  to  detector  A.   Phase  sensitive 
detection  of  the  first  or  third  harmonic  would  yield  an 
error  signal.   The  second  harmonic  is  unsuitable  for  error 
signal  production. 

F.   THE  SERVO  LOOP  INTEGRATOR 

It  is  desired  to  maintain  the  P(20)  line  in  coincidence 
with  the  center  frequency  of  a  Stark  shifted  NH?D  transition 
as  accurately  as  possible.   When  the  two  transitions  are  in 
coincidence  the  error  signal  is  zero.   Using  Fig.  19  as 
an  example,  consider  a  slight  temperature  increase  in  the 
cooling  water.   The  laser  cavity  will  expand  and  pull  the 
oscillation  to  a  lower  frequency  according  to  (45) .   This 
results  in  a  negative  o  and  a  positive  error  signal  pro- 
portional to  <5   (assuming  6  >  -.2  in  Fig.  19).   Feeding  the 
error  signal  directly  to  the  PZ  supply  then  causes  the 
supply  to  increase  the  PZ  voltage  and  make  a  cavity  length 
contraction  correction.   But  as  the  cavity  length  returns 
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to  the  original  length  the  error  signal  decreases  until  an 
equilibrium  point  is  reached  at  some  finite  nonzero  <5 ,  such 
that  the  frequency  difference  6  produced  by  the  PZ  voltage 
produces  just  the  right  error  signal  to  produce  that  PZ 
voltage.   It  is  obvious  that  6  could  never  return  to  zero 
as  the  error  signal  would  be  zero  indicating  no  length 
correction. 

A  method  is  needed  to  maintain  the  required  PZ  voltage 
as  the  error  signal  returns  to  zero.   This  could  be  accom- 
plished by  integrating  in  time  the  error  signal  produced 
by  the  lock-in  and  producing  a  new  error  signal  to  drive 
the  PZ  supply  which  is  proportional  to  the  total  time 
integral  of  the  first  error  signal.   To  the  author's  know- 
ledge, the  use  of  an  integrator  in  the  servo  loop  has  not 
been  mentioned  before  in  the  literature. 

The  integrator  is  placed  between  the  lock-in  amplifier 
and  the  PZ  power  supply  as  indicated  in  Fig.  6.   Now  as 
the  error  signal  begins  to  return  toward  zero,  the  integrator 
continues  to  increase  its  output  until  the  error  signal  is 
actually  zero.   At  5  =  0  the  integrator  maintains  a 
constant  D.C.  output  indefinitely  or  until  a  nonzero  error 
signal  indicates  a  need  to  change.   If  the  error  signal 
again  went  positive  the  integrator  output  increases  still 
further,  until  a  return  to  zero  of  the  error  signal  indi- 
cates for  the  integrator  to  maintain  a  constant  output  at 
the  new  voltage  reached. 


If  the  error  signal  went  negative  the  integrator  would 
lower  its  output  and  continue  to  do  so  until  the  error 
signal  actually  returned  to  zero. 

The  ideal  integrator  would  respond  to  both  rapid  and 
slow  changes  of  the  error  signal.   It  would  integrate 
rapidly  enough  so  as  to  stabilize  against  acoustic  vibrations 
yet  also  be  able  to  accommodate  slow  changes  in  tempera- 
ture.  It  should  be  designed  to  stabilize  over  the  maximum 
temperature  change  expected  and  maximum  pressure  change 
expected.   Such  a  device  would  produce  a  highly  stable 
system  suitable  for  operation  in  the  field  or  in  the 
laboratory. 

G.   THE  STARK  CELL 

Stabilization  was  never  achieved  due  to  numerous  design 
problems  in  the  Stark  cell.  Six  different  cells  were  made. 
These  are  discussed  below. 

1.   Stark  Cell  One 

The  first  Stark  cell  had  stainless  steel  plates 
10  cm  long  and  1 . 5  cm  wide  glued  with  plastic  spacers  into 
a  plastic  tube  1.7  cm  I.D.   Stainless  leads  fastened  to 
the  plates  were  sealed  with  glue  through  holes  in  the  side 
of  the  tube.   The  cell  was  made  to  accommodate  one- inch 
NaCl  windows  as  were  all  the  cells.   The  plate  material  was 
too  thin  and  as  the  glue  hardened  on  the  leads,  the  plates 
warped. 
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2.  Stark  Cell  Two 

This  was  made  with  heavy  stainless  plate  10  cm 
long  and  3  cm  wide  with  a  plate  separation  of  1.4  mm. 
The  maximum  observed  D.C.  modulation  was  20%  at  a  D.C. 
voltage  of  490  volts  and  a  pressure  of  4.7  torr.   The 
electric  field  was  perpendicular  to  the  laser  polarization. 
The  problem  was  that  breakdown  occurred  at  500  volts. 
Also  diffraction  was  a  problem  with  the  narrow  plate 
separation. 

3.  Stark  Cell  Three 

This  was  made  with  plates  20  cm  long,  3 . 5  cm  wide 
and  about  3  mm  thick.   Care  was  taken  to  round  the  sides 
and  corners  of  the  plates.   The  plates  were  separated  to 
4  mm  in  an  effort  to  eliminate  breakdown  and  minimize 
diffraction.   The  cell  had  a  gross  air  leak  which  could 
only  be  remedied  by  redesign.   The  cause  of  the  leak  was 
the  rubber  "0"  rings. 

4.  Stark  Cell  Four 

This  was  the  same  as  above  with  slightly  modified 
end  window  mounts.   Slight  modulation  was  obtained  at  1400 
volts  at  a  pressure  of  47  torr  with  the  E  field  perpendicular 
to  P.   The  intensity  through  the  cell  was  .12  W.   Breakdown 
also  occurred  at  this  potential.   At  lower  pressures 
breakdown  occurred  at  lower  potentials . 

It  was  then  pointed  out  that  in  a  low  pressure  gas 
electricity  does  not  always  take  the  shortest  path.   In 
fact,  it  was  observed  in  this  cell  that  at  certain  pressures 
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the  discharge  occurred  between  the  outside  two  surfaces  of 
the  plates  where  the  distance  was  greater  rather  than 
between  the  two  inside  surfaces.   It  is  a  fact  that  break- 
down potential  is  not  a  function  of  distance  but  of 
pressure  times  distance.   A  plot  of  breakdown  potential  as 
a  function  of  p  x  d  is  a  skewed  parabola  with  a  minimum 
at  some  value  of  p  x  d.   This  value  is  called  the  pd 
minimum  or  Paschen  minimum.   To  the  right  of  the  pd  minimum 
breakdown  potential  increases  with  increasing  pressure  and 
distance  as  would  normally  be  expected.   But  to  the  left 
of  the  pd  minimum  the  breakdown  potential  decreases  with 
increasing  pressure  and  distance.   Consequently  if  one  is 
unfortunate  in  his  choice  of  pressure  and  distances  such 
that  he  is  operating  to  the  left  of  the  pd  minimum,  increasing 
the  plate  separation  will  only  make  matters  worse. 
5.   Stark  Cell  Five 

The  above  cell  was  then  disassembled  and  mylar  was 
glued  to  the  inside  surfaces  of  the  plates.   All  other 
exposed  metal  surfaces  including  the  lead  wires  were 
coated  with  RTV  silicone  rubber.   As  the  cell  was  placed 
under  vacuum  the  glue  and  RTV  began  to  outgas  forming  gas 
pockets  on  the  plates  and  in  the  RTV.   The  required  1500 
volt  potential  placed  on  the  plates  punched  through  the 
RTV  and  still  resulted  in  breakdown. 

The  problem  was  reevaluated.   The  operating  pressure 
could  not  be  changed  enough  to  solve  the  breakdown  problem 
as  the  Stark  effect  ceases  to  give  good  modulation  at  too 
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high  or  too  low  a  pressure  [Ref.  18].   Two  physical 
characteristics  were  desired  which  appeared  to  be  mutually 
exclusive:   a  wide  channel  between  the  plates  for  ease  of 
alignment  and  a  minimum  of  diffraction,  and  the  elimination 
of  the  breakdown  problem. 

6.   Stark  Cell  Six 

A  completely  new  design  approach  was  taken:   to 
effectively  place  all  metal  surfaces  outside  of  the  cell, 
isolated  from  the  low  pressure  gas  filled  channel.   A 
3  mm  by  10  mm  wide  channel  was  cut  in  a  solid  22  cm  long 
block  of  plastic.   The  plates  were  embedded  in  the  plastic 
on  either  side  of  the  3  mm  width.   The  design  called  for 
1  mm  of  plastic  between  each  plate  and  the  channel,  resulting 
in  a  plate  separation  of  5  mm.   Due  to  crazing  of  the  1  mm 
thick  wall  in  the  machining  process,  a  film  of  mylar  was 
added  to  the  underside  of  each  plate  when  they  were  embedded 
in  the  plastic.   This  resulted  in  an  uneven  plate  separation 
of  approximately  7  mm.   The  cell  was  tested  to  3000  volts 
with  no  breakdown.   It  also  held  a  good  vacuum.   However, 
it  is  anticipated  that  the  uneven  plate  separation  may 
result  in  a  broadened  or  distorted  absorption  profile. 
Testing  was  not  carried  out  further  due  to  a  permanent 
failure  of  the  laser  and  its  power  supply. 

This  design  appears  to  be  the  best  approach.   With  this 
plate  separation  of  7  mm,  wider  plates  should  be  used  to 
minimize  fringe  effects  of  the  electric  field.   Also  the 
fabrication  process  needs  to  be  modified  so  that  the  plates 
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are  parallel,  but  sealed  in  plastic  with  no  crazing  in  the 
high  field  areas.   The  plates  should  be  thick  so  that  the 
corners  and  edges  can  be  rounded  to  minimize  possible 
breakdown  of  the  plastic.   One  solution  might  be  to  first 
cement  a  heavy  plate  to  a  block  of  plastic  then  machine  the 
channel  and  cement  the  two  halves  together. 
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IV.   CONCLUSION 

The  instabilities  of  the  Sylvania  C02  Model  941  laser 
have  been  studied  and  the  laser  determined  to  be  suitable 
for  Stark  cell  stabilization.   It  was  found  that  the  exist- 
ing unregulated,  unfiltered  laser  power  supply  was  unsuitable 
for  a  stabilized  system  and  a  new  well  regulated  power 
supply  needs  to  be  purchased.   Laser  characteristics 
required  to  build  a  stabilized  system  such  as  piezoelectric 
parameter  and  thermal  expansion  coefficient  have  been 
measured.   The  optimum  discharge  current  was  found  to  be 
7.5  mA.   All  the  physical  and  operating  characteristics  of 
the  laser  pertinent  to  the  design  of  a  stabilized  system 
have  been  studied  and  listed. 

A  computer  analysis  has  been  carried  out  to  demonstrate 
the  feasibility  and  production  of  an  error  signal  with 
existing  equipment.   Both  a  Pyroelectric  and  HgCdTe  detector 
are  suitable  for  error  signal  production  by  phase  sensitive 
detection  of  the  first  harmonic.   A  secondary  result  of  the 
computer  analysis  was  the  demonstration  of  the  need  for  a 
servo  loop  integrator.   The  purpose  and  operation  of  such 
an  integrator  was  discussed.   Theoretically  the  integrator 
should  provide  increased  stability  over  that  indicated 
thus  far  in  the  literature. 

A  Stark  cell  with  no  known  major  design  flaws  has  been 
constructed  and  tested  to  3000  volts  with  no  breakdown.   It 
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is  vacuum  tight  to  at  least  10~   torr.   It  is  expected  to 
be  successfully  used  in  the  stabilization  system. 

All  preliminary  engineering  studies  have  been  carried 
out  for  assembling  a  Stark  stabilized  C02  laser  system 
operating  on  the  P(20)  mode.   This  system  will  be  opera- 
tional in  both  the  laboratory  and  the  field. 
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APPENDIX  A 
MODEL  941  SYLVANIA  LASER  SPECIFICATIONS 

Output  wavelength:   10.6  ym 

Power  output :   3  watts 

Beam  divergence:   4  mrad 

Beam  diameter  to  1/e  point:   4  mm 

Beam  polarization:   vertical 

Cavity  length:   45  cm 

Bore  diameter:   ~  6  mm 

Cooling  requirements:   tap  water  @  15  to  30  ml/s 

Output  mirror:   70%  reflecting,  dielectric  flat 

Piezoelectric  mirror:   100%  reflecting,  300  cm  radius  of 

curvature 

Piezoelectric  ^  parameter 

As  given  by  Sylvania:   -4  x  10   urn/volt 

Measured:   -3.8  x  10"   ym/volt 
Laser  voltage:   14  kV  @  7.6  ma 
Tube  voltage:   9.5  kV  §  7.6  ma 
Series  ballast  resistance:   600  kQ, 
Equivalent  tube  resistance  §  7.6  ma:   1.2  7  MA 
Pressure  broadened  linewidth:   100  MHz 
Doppler  broadened  linewidth:   50  MHz 

ij-  mode  separation:   333  MHz 
Cavity  linewidth:   19  MHz 
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APPENDIX  B 
STARK  CELL  SPECIFICATIONS 

Two  parallel  plates  totally  encased  in  plastic  with  a 
3  mm  x  10  mm  light  beam  channel  at  the  axis  of  the  cell 
between  the  plates.   Channel  is  sealed  off  with  one  inch 
NaCl  windows  at  each  end  and  is  vacuum  tight.   Channel  is 
filled  through  one  port. 

Plate  separation:   7  mm 
Plate  width:   3.6  cm 
Plate  length:   20.3  cm 
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APPENDIX  C 
HP  9810A  CALCULATOR  PROGRAM 
FOR  THE  VOIGT  PROFILE 

The  Voigt  profile  is  the  convolution  of  a  Lorentzian 
and  a  Gaussian  profile.   This  HP  9810A  program  convolves 
numerically  a  Gaussian  g-,  (x)  and  a  Lorentzian  g2Cx)  which 
have  equal  heights  at  the  center  frequency,  i.e., 


MV  =  MV  •  <C1) 


Solving  (CI)  and  setting  vQ  =  0  and  j^-   =  -  the  Gaussian 
and  Lorentzian  profiles  are  simplified  to 


2 
gx(x)  =  exp  [-ttx  ] 

and  (C2) 

g  (x)  =  T- 

2  [1  +  (HX)Z] 


These  were  the  equations  plotted  in  Fig.  5.    Equations  (C2) 
are  convolved  according  to 


x 
max 

§(X)  V   gl(7)  S2(X  '  7)  d7  (C3) 

min 
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The  flow  diagram  and  complete  program  are  given  on  the  follow- 
ing pages.   The  counter  containing  the  value  of  n  is 
register  (000)  .   The  value  of  y  is  stored  in(b),  the  value 
of  x  in  (a).   dx  and  dy  are  stored  in  (100).   1/tt   is  stored 
in  (099).   (098)  is  a  temporary  storage  for  holding  the 
previous  value  of  f  for  the  numerical  integration  routine. 
Registers  (001)  through  (093)  are  reserved  for  storing  S 

where  S   goes  into  (n) .   x    is  stored  in  (097) . 
n  max 

x    -  x  . 

=  ma* EiH  is  stored  in  (096).   n  +  1  is  stored 

max       dx 

in  (095)  and  A  which  is  the  normalizing  factor  is  stored 

in  (094).   All  registers  should  be  cleared  before  beginning. 

Begin  by  prestoring  these  typical  values: 

dx  =  dy  =  .1  ■»■  (100) 

x    =  -x  .   =  4  -»■  (097) 
max     mm 

n    =  80  -  (096) 
max 

y  =  -4  ->  (b) 
x  ■  -4  ■+■  (a) 

(^.)2  =  "ij  -  (099) 

IT 

Set  up  the  plotter  for  xmax  =  4,  *min  =  -4  and  ymax  =  1, 

y    =0   Load  the  program  and  press  continue.   The  computer 
7min 

will  plot  out  the  Voigt  profile  in  approximately  20  minutes. 
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APPENDIX  D 
HP  9810A  CALCULATOR  PROGRAM  FOR 
STARK  FREQUENCY  DISCRIMINANT 

I.   FREQUENCY  DISCRIMINANT  FOR  PYROELECTRIC  DETECTOR 

The  following  program  plots  b   in  Eq.  (66)  as  a  function 
of  5,  where  f(t)  is  given  by  (64)  with  S  =  1  and  Av  =  1 . 
The  first  step  in  the  program,  K2,  sets  the  calculator  to 
work  in  radians.   5  is  prestored  in  (000),  Av f  is  prestored 
in  (006)  ,  and  n  (which  determines  the  harmonic)  is  pre- 
stored in  (007).   The  values  of  the  constants  A,  B,  C,  D, 
and  E  in  (64)  are  calculated  and  stored  in  registers  (001)  , 
(002),  (003),  (004),  and  (005),  respectively,  in  the 
program.   The  variable  t  is  stored  in  (a)  and  dt  in  (b) . 
(008)  is  a  temporary  store  for  the  numerator  in  (64) . 
(012)  is  a  temporary  store  for  the  numerical  integration 
routine.   (010)  accumulates  the  value  of  the  integral  in 
(66). 

Begin  with  all  registers  cleared  and  the  program  loaded. 
Set  up  plotter  for  -1  <  x  <  1  and  -2  <_  y  <_  2  or  some  other 
appropriate  limits.   Prestore  -1  in  register  (000),  .5  in 
(006),  1  in  (007),  .01  or  .05  in  (b) .   Press  continue  and 
Fig.  20  will  be  plotted  out  from  left  to  right. 
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II.   FREQUENCY  DISCRIMINANT  FOR  HgCdTe  DETECTOR 

The  following  program  plots  a   in  Eq.  (65)  as  a  function 
of  5,  where  f(t)  is  given  by  (61)  with  S  =  1  and  Av  =  1.   As 
before  (000)  ,  (006)  ,  (007)  ,  and  (b)  contain  the  appropriate 
prestored  values.   In  the  program  (011)  is  the  numerical 
integration  routine  temporary  storage  register  and  (009) 
accumulates  the  value  of  the  integral. 
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APPENDIX  E 
LOW  PASS  L-SECTION  FILTER 

Below  is  a  schematic  of  the  simple  filter  added  to  the 
high  voltage  D.C.  laser  power  supply  to  reduce  the  60  Hz 
current  ripple.   The  filter  consisted  of  a  .5  microfarad 
16  kilovolt  capacitor  and  a  45  Henry  choke  which  was 
actually  the  secondary  winding  of  a  13  kV  transformer. 
Since  the  input  voltage  was  15  kV,  the  transformer,  being 
insulated  for  2000  volts  less,  was  isolated  from  ground  on 
plastic  supports.   The  filter  successfully  reduced  the  60  Hz 
ripple  in  the  power  supply  output  by  a  factor  of  50. 


"15  KV  DC 


45  H 

UNFILTERED 
10%    RIPPLE 


.5^7- 


TO    LASER 
-o    > 


.2%    RIPPLE 
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APPENDIX  F 
COPPER  CONSTANTAN  THERMOCOUPLE  TEMPERATURE  VS.  E.M.F 

The  following  table  gives  values  of  thermocouple 
temperature  and  voltage  with  reference  to  ice  water. 

E.M.F.  (mV)  Temperature  (°C) 


0 

0 

.1 

2.61 

.2 

5.17 

.3 

7.72 

.4 

10.28 

.5 

12.78 

.6 

15.33 

.7 

17.83 

.8 

20.33 

.9 

22.78 

1.0 

25.28 

1.1 

27.72 

1.2 

30.11 

1.3 

32.56 

1.4 

35.00 

1.5 

37.39 

1.6 

39.78 

1.7 

42.11 

1.8 

44.50 

1.9 

46.83 

2.0 

49.17 

109 


APPENDIX  G 
RELATIONSHIP  OF  TUBE  WALL  TEMPERATURE 
TO  COOLING  WATER  OUT  TEMPERATURE 

The  laser  discharge  region  between  the  cavity  mirrors 
is  a  pyrex  glass  tube  .7  cm  O.D.  and  34.3  cm  long  sur- 
rounded by  a  glass  water  jacket  3  cm  I.D.   It  is  desired 
to  find  the  average  wall  temperature  of  the  .7  cm  discharge 
tube  knowing  the  bulk  water  temperature  out.   Reference  24 
is  used  for  this  calculation.   The  problem  is  an  annular 
tube  flow  type.   It  is  assumed  that  the  outer  water  jacket 
tube  is  at  the  bulk  water  temperature  and  that  no  heat 
flows  through  it  to  the  surrounding  air. 

The  following  parameters  are  given,  through  experimental 
measurement  or  a  handbook  of  physical  constants. 

flow  rate:   14.6  ml/s 

inlet  temp . :   22. 1°C 

outlet  temp.:   23.8°C 

voltage  across  tube  and  balast  resistor:   14  kV 

voltage  across  tube:   9 . 5  kV 

tube  current:   7.6  lA 

laser  power  output:   3  watts 

specific  heat  of  water:   Cp  =  1  cal/g-C° 

/   3 
density  of  water:   p  =  1  g/cm 

viscosity  of  water  at  22°C:  \i   =  .009403  g/cm-s 
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thermal  conductivity  of  water  at  0°C:   k   =  .001348  cal/cm-s-°C 


w 


thermal  conductivity  of  water  at  40°C:   k  =  .001499  cal/cm-s-°C 


w 


thermal  conductivity  of  pyrex:   kg  =  .0026  cal/cm-s-°C 

hydraulic  diameter  2(r   -  r.):   D  =  2 . 3  cm 

o    1/    e 

tube  surface  area"  Ag  =  tt  ( .  7  cm)  (34.3  cm)  =75.4  cm2 


1.  Determine  if  the  flow  is  laminar  or  turbulent  by  calcu- 
lating the  Reynolds  number. 

2      2  ? 

cross  sectional  area  of  water:   A  =  7T(r   -  r.  )  =  6.68  cm 

o      1 

velocity  of  water:   V  =  14,6  ml^s  =2.19  cm/s 

6.68  cnT 

Re  =  p!^  =  U  g/cmll(oh\9   ™£sULl   cm)  =  536 
u  .009403  g/cm-s 

.  .  Flow  is  definitely  laminar. 

2.  Determine  Nusselt  number  for  laminar  flow  in  a  circular 
tube  annulus,  using  [24]. 


r  . 
r 


r*  =  -i  =  .23 


o 


Nu-  =  8.212 


3.   Find  the  heat  added  to  the  water  by  the  laser  discharge 
tube  only.   This  is  approximately  equal  to  the  voltage 
across  the  tube  times  the  current  minus  the  power  output. 

SJS  =(9.5  x  103)(7.6  x  10"3)  -  3  =  69.2  watts  =  16.54  cal/s 
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4.   Find  the  temperature  difference  AT  across  the  tube  and 
the  bulk  water  temperature  given  the  temperature  out. 

dQ  =  r      at  ^ 
oT   Cp  T  oT 


16.54  cal/s  =  (1  |ijU-)  AT  (14.6  g/s) 


AT  =  1.1  C' 


T   +  T.    T   +  (T   -  1.1)        ,  , 
T   _  _o l    o    ^  o '    T    1 . 1 

b      2    '         2       '   o  "  ~T~ 


5.   The  surface  coefficient  of  the  inner  tube  is  found 
from  the  Nusselt  number  according  to 


k  Nu. 

u   _   w    1 


6.   The  rate  of  heat  transfer  through  a  . 5  mm  glass  tube 
wall  with  an  inside  surface  temperature  Tw.  and  an  out 
side  surface  temperature  Tw   is  given  by 


jn  Tw.  -  Tw 

dQ  =  k  a  r  1         si 

Jt        KgAs  l  .05  cm  J 


16.54  S|i  =  (.0026  _^|^_)(75.4  cm2)  (Twi  -  Tws)/.05  cm 


Tw.  -  Tw   =  4.2  C 
l     s 
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7.   Derive  a  general  formula  for  the  average  tube  wall 

temperature  given  that  the  rate  of  heat  transfer  through 
the  tube  wall  equals  the  rate  through  the  surface 
layer  equals  the  rate  dissipated  by  the  water. 


An  Tw.     -    Tw  , 

§  -  kg  As     '.os      -  hs  As  (T»s  -  V  =  S  AT  a* 


Tw.    +    Tw 
Tw 


A       2 
From  the  above  it  is  easily  shown  that 

T     1  r.05  .  dQ    ,  1   .  dQ   T 

TwA  =  2  (k-A7}  dT  +  (h7T)  dT  +  Tb 

Substitution  of  the  equations  for  T,  and  hs  in  4  and  5 
into  this  equation  yields 


D 


Tw  1    f  05    ,    dQ         f      De         ,    dQ        T 

TwA   "   I   CI^AT:)    3t        Ck   Mu-Aj    dT  o 


1.1 


or 


w 


=  45.8°  +  TQ  @  0°C 


=  41.4°  +  T   @  40°C 
Clearly  this  proves  that  ATw.  =  AT 
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APPENDIX  H 
PHOTOGRAPHS  OF  EQUIPMENT 

Figure  (H-l)  is  a  photograph  of  the  whole  experimental 
setup.   The  L-section  filter  is  on  top  of  the  laser  power 
supply  in  the  lower  right-hand  corner.   The  cooling 
water  temperature  control,  filter,  and  regulator  are  at 
the  sink  in  the  upper  right-hand  corner. 

Figure  (H-2)  is  a  photograph  of  the  oscilloscope  and 
rack-mounted  instruments.   From  top  to  bottom  the  instru- 
ments in  the  rack  are:   audio  power  amplifier,  lock-in 
amplifier,  3kV  D.C.  power  supply,  2kV  D.C.  power  supply, 
audio  oscillator,  and  nonovoltmeter . 

Figure  (H-3)  is  a  photograph  of  the  laser  and  optical 
components.   The  laser  is  middle  right  in  the  photograph. 
Note  the  thermometer  and  thermocouple  in  the  coolant  tube 
coming  out  of  the  rear  of  the  laser.   The  gas  bottles  and 
dewar  for  the  thermocouple  ice  bath  are  at  the  bottom  of 
the  photograph.   From  the  laser  counterclockwise  are:   the 
chopper,  beam  splitter,  spectrum  analyzer,  Stark  cell,  and 
pyroelectric  detector. 

Figure  (H-4)  is  a  photograph  of  the  Stark  cell.   From 
left  to  right  are:   the  beam  splitter,  chopper  motor,  laser, 
Stark  cell,  and  detector.   Note  that  the  chopper  motor  is 
acoustically  isolated  from  the  table  with  a  metal  bar  sup- 
ported directly  to  the  floor. 
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FIGURE     H-2 
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FIGURE    H-3 
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